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EXECUTIVE  SUMMARY 


A.  OBJECTIVE 

The  objective  of  this  guidance  manual  is  to  provide  a  tool  to  U.S.  Air  Force  engineers  and 
decision  makers  for  use  in  evaluation  and  selection  of  air  stripping  with  emissions  control  technology. 

B.  BACKGROUND 

The  development  of  this  manual  was  part  of  a  larger  activity  for  the  U.S.  Air  Force  focused  on 
the  field  study  and  demonstration  of  the  state  of  the  art  in  air  stripping  with  emissions  control 
technology.  The  intent  of  the  guidance  manual  is  to  provide  managers  and  engineers  with  necessary 
information  so  that  decisions  relating  to  the  application  of  this  technology  can  proceed  on  a  rational 
basis.  Although  the  information  provided  in  the  manual  is  available  from  a  number  of  sources,  it  had 
not  been  organized  so  that  it  could  be  readily  used  for  air  stripping  technology  evaluation,  especially 
by  one  not  very  skilled  in  the  art.  In  this  guidance  manual,  the  information  is  organized  and  provided 
to  the  user  in  an  interactive  spreadsheet  format. 

C.  SCOPE 

The  results  of  the  field  tests  for  the  U.S.  Air  Force,  along  with  other  relevant  information,  were 
translated  into  a  computer  model  used  to  assess  the  design  and  economic  aspects  of  relevant 
technology  and  operating  parameters  for  groundwater  remediation  by  air  stripping  with  emissions 
control.  The  technology  basis  for  the  computer  model  and  operating  instructions  for  the  spreadsheet 
are  included  in  the  guidance  manual.  Also,  the  spreadsheet  was  used  to  evaluate  example  cases  of 
air  stripping  with  emissions  control.  The  results  of  the  calculations  are  presented  and  explained  in 
the  manual. 

D.  METHODOLOGY 

The  approach  used  in  this  activity  was  to  provide  a  useful  computer  model  in  an  easy  to  use 
format.  The  economic  simulation  was  developed  as  spreadsheets  that  may  be  used  on  a  personal 
computer.  Multiplan  by  Microsoft  was  the  software  chosen  because  of  its  capability  to  automatically 
perform  iterative  calculations  (needed  when  circular  dependencies  exist  between  spreadsheet  cells). 
In  the  development  of  the  computer  model,  information  was  utilized  from  the  literature  and  also  from 
the  results  of  the  field  study  of  air  stripping  with  emissions  control  conducted  for  the  U.S.  Air  Force. 
E  TEST  DESCRIPTION 

The  Eglin  Air  Force  Base  air  stripping  field  study  is  described  in  full  in  ESL  TR  90-51,  Air 
Stripping  and  Emissions  Control  Technologies:  Field  Testing  of  Countercurrent  Packings.  Rotary  Air 

Hi 


Stripping.  Catalytic  Oxidation,  and  Adsorption  Materials.  Results  from  this  study  were  used  in  the 
development  of  the  computer  model  for  this  guidance  manual. 

F.  RESULTS 

Analysis  of  the  lifetime  operating  costs  for  application  of  air  stripping  technology  for  remediation 
of  contaminated  groundwater  indicates  that  the  use  of  emissions  control  devices  for  systems  for  air 
stripping  of  VOCs  from  groundwater  considerably  increases  the  costs  of  such  operations,  cost 
considerations  favor  the  use  of  lower  values  of  the  stripping  factor,  and  the  costs  of  the  remediation 
of  contaminated  groundwater  are  largely  controlled  by  operating  rather  than  capital  costs. 
a  CONCLUSIONS 

The  spreadsheet-based  model  developed  and  described  herein  can  be  readily  utilized  by  Air 
Force  personnel  for  the  evaluation  and  selection  of  an  air  stripping  with  emissions  control  technology 
for  the  remediation  of  VOC  contaminated  groundwater.  If  desired,  the  correlations  for  design  and 
cost  estimations  can  be  updated  as  new  information  becomes  available.  By  varying  financial  input 
parameters,  such  as  depreciation  period,  interest  rates,  etc.,  different  scenarios  may  be  evaluated  for 
particular  operating  conditions  of  groundwater  flowrate,  contaminant  concentration,  etc. 
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PREFACE 


This  report  was  prepared  by  the  Oak  Ridge  National  Laboratory,  P.  0.  Box 
2008,  Oak  Ridge,  TN  37831-6044,  for  the  Air  Force  Engineering  and  Services 
Center,  Engineering  and  Services  Laboratory,  Tyndall  Air  Force  Base,  Florida,  as 
partial  fulfillment  of  the  statement  of  work  entitled  "Air  Stripping  with 
Emissions  Control"  in  accordance  with  DOE  Interagency  Agreement  No.  1489-1489-Al . 
The  period  of  performance  of  the  work  was  from  April  1987  to  June  1989.  Oak 
Ridge  National  Laboratory  is  managed  by  Martin  Marietta  Energy  Systems,  Inc.,  for 
the  U.S.  Department  of  Energy  under  contract  DE-AC05-840R21400. 

This  document  details  the  results  of  activities  performed  under  Task  4.4  of 
the  statement  of  work.  Phase  4:  IRP  Pha.se  IV  Users  Manual.  A  related  document 
completed  under  the  same  contract  is  ESL  TR  90-51,  Air  Strippino  and  Emissions 
Control  Technologies:  Field  Testing  of  Countercurrent  Packings.  Rotary  Air 
Stripping,  Catalytic  Oxidation,  and  Adsorption  Materials.  The  AFESC/RDVW  Project 
Officers  for  this  effort  were  Captain  R.  A.  Ashworth  and  Captain  Edward  Marchand. 

Mention  of  trade  names  or  commercial  products  within  this  document  does  not 
constitute  endorsement  or  recommendation  for  use. 

Copies  of  the  V0C-2.100  and  V0C-2.300  spreadsheets  are  available  from  the 
Chemical/Physical  Treatment  Technology  Area  Manager,  HQ  AFESC/RDVW,  Tyndall  AFB, 
FL.  Include  a  diskette  with  your  request. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be 
available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

EDWARD  G.  MARCHAND,  Capt,  USAF,  BSC  FRANK  P.  GALLA^ER  III,  Col,'  USAF 

Chemical/Physical  Treatment  Director,  AF  c4vil  Engineering 

Technology  Area  Manager  Laboratory 


Chief,  Environics  Division 
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MANUAL  FOR  ESTIMATING  COST  OF  VOC  REMOVAL  FROM 
GROUNDWATER  CX)NTAMINATED  WITH  JET  FUEL 

SECTION  I 
INTRODUCTION 


A.  OBJECTIVE 

This  manual  is  intended  to  provide  relevant  information  for  persons  involved  in  the  remediation 
of  groundwater  contaminated  with  jet  fuel  or  similar  compounds.  Developing  strategies  for 
groundwater  remediation  consists  of  four  principal  steps:  (1)  site  characterization;  (2)  determination 
of  treatment  requirements;  (3)  remediation  system  design;  and  (4)  system  fabrication,  installation,  and 
operation.  This  manual  is  designed  to  be  of  primary  benefit  for  item  (3)  of  the  above  list;  it  may  also 
be  of  benefit  for  items  (2)  and  (4).  Typical  remediation  technologies  are  identified  and  discussed; 
the  costs  of  installed  systems  based  on  these  technologies  may  be  estimated  using  the  enclosed 
software  and  an  IBM-compatible  personal  computer.  Much  of  the  information  provided  in  this 
manual  is  from  a  field  study  of  air  stripping  with  emissions  control  performed  by  the  Oak  Ridge 
National  Laboratory  (ORNL)  for  the  United  States  Air  Force  (Reference  1). 

Jet  fuel  contains  a  large  percentage  of  soluble  hydrocarbon  compounds  that  are  easily  removed 
from  groundwater  by  aeration;  these  compounds  are  called  volatile  organic  compounds  (VOCs)  and 
have  a  Henry’s  coefficient  of  greater  than  16  to  32  ft’  atm/pound  mole  [1  to  2  m’  atm/kmol]. 
Benzene  is  a  VOC,  having  a  Henry’s  coefficient  estimated  to  be  76  ft’  atm/pound  mole  [4.7  m’ 
atm/kmol]  at  68®F  (Reference  2);  naphthalene  is  not  considered  to  be  a  VOC,  having  a  Henry’s 
coefficient  estimated  to  be  6.7  ft’  atm/pound  mole  (0.42  m’  atm/kmol]  at  68  “F  (Reference  3). 
Compounds  with  Henry’s  coefficient  of  less  than  16  ft’  atm/pound  mole  [1  m’  atm/kmol],  such  as 
naphthalene,  are  not  ordinarily  removed  very  efficiently  in  aeration  devices  designed  for  VOC 
removal. 

Stripping  fuel  components  from  water  merely  transfers  them  to  another  medium,  usually  air. 
Abatement  of  gaseous  hydrocarbons  by  adsorption  onto  activated  carbon  has  been  utilized  for  many 
years;  this  process  transfers  the  contaminants  to  another  medium  —  activated  carbon  —  which 
requires  disposal  or  regeneration.  Abatement  by  catalytic  oxidation  potentially  offers  a  very  high 
destruction  of  fuel  components,  producing  COj  and  H2O. 
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Infonnation  contained  in  this  manual  and  accompanying  software  provides  the  user  with  a 
method  of  estimating  unit  processing  costs  (U.S.  dollars/1000  gallons)  for  remediation  of  groundwater 
contaminated  with  jet  fuel  components.  Two  options  are  offered  for  air  stripping  —  a  traditional 
packed  tower  and  a  rotary  air  stripper.  Also,  two  options  are  offered  for  emissions  control  — 
activated  carbon  adsorption  and  catalytic  oxidation.  Estimated  treatment  costs  include  fixed  capital 
investment  recovery,  direct  and  indirect  operating  costs,  and  a  method  of  accounting  for  uncertainty 
in  any  cost  component  The  estimates  in  these  primary  cost  categories  are  used  as  inputs  to  an 
operating  lifetime  cost  analysis.  The  cleanup  costs  are  usually  dominated  by  annual  operating  costs 
rather  than  fixed  capital  investments.  Cost  estimating  techniques  are  compatible  with  procedures  of 
the  U.S.  General  Accounting  Office.  American  units  (pounds-feet-seconds)  are  used  in  this  report. 
When  not  typical  of  supporting  information,  other  more  appropriate  units,  usually  metric,  are  given 
in  brackets,  [  ]. 

B.  BACKGROUND 

Contamination  of  groundwater  by  VOCs  is  an  important  environmental  problem.  Sources  of 
such  contamination  include  accidental  spills  and  leaking  storage  tanks  and  transport  lines. 
Contamination  of  soil  and  groundwater  from  a  leaking  storage  tank  is  illustrated  in  Figure  1.  As 
shown  in  this  figure,  remediation  of  a  contaminated  aquifer  is  part  of  a  larger  problem  of  remediating 
the  contaminated  site,  which  includes  the  unsaturated  zones.  Site  remediation  usually  involves 
containment  of  the  contaminant(s).  The  National  Oil  and  Hazardous  Substances  Contingency  Plan 
(Reference  4)  identifies  three  general  categories  of  remediation  activities:  (1)  initial,  (2)  source 
control,  and  (3)  off-site  measures.  Pumping  of  groundwater  can  be  useful  for  all  three  categories. 
Groundwater  that  ordinarily  flows  continuously  through  the  site  is  an  important  agent  for  transporting 
soluble  contaminants.  Pumping  activities  usually  involve  depressing  the  groundwater  table  in  the 
affected  region  by  pumping  from  a  system  of  wells.  Contaminated  water  from  such  an  operation 
requires  treatment,  as  necessary,  for  its  discharge;  such  treatment  is  the  subject  of  this  manual.  A 
block  diagram  for  the  remediation  of  contaminated  groundwater  by  air  stripping  with  emissions 
control  is  presented  in  Figure  2. 

A  review  of  relevant  technology  was  recently  presented  by  Singh  and  Counce  (Reference  5). 
Air-water  contactors  may  be  designed  to  give  high  removal  of  VOCs  from  water;  contaminants  of  less 
volatility  will  be  removed  simultaneously  to  a  lesser  extent.  The  contaminated  effluent  air  from  the 
air-water  contactor  is  often  a  matter  of  concern.  If  this  stream  is  not  suitable  for  discharge,  then 
treatment  to  restore  acceptable  air  quality  will  be  required.  During  the  remediation  activity,  the 
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Figure.  1.  Subsurface  Behavior  of  Spilled  Hydrocarbons. 
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concentration  of  contaminants  in  the  groundwater  will  be  reduced  to  an  acceptable  level  for 
discharge.  At  some  point  in  this  activity,  the  effluent  air  quality  may  improve  to  the  point  that  any 
emissions  control  originally  required  is  no  longer  necessary.  The  objective  of  an  air  stripper  with  an 
emissions  control  unit  is  to  remove  VOC  contaminants  from  water  so  that  the  water  is  suitable  for 
discharge  for  the  life  of  the  remediation  activity;  a  corresponding  objective  is  to  operate  the  unit  so 
that  all  other  appropriate  environmental  regulations  are  observed  for  the  same  time  period. 

C  SCX)PE/APPROACH 

The  purpose  of  this  manual  is  to  provide  a  means  of  assessing  the  unit  processing  costs  (U.S. 
dollars/1000  gallons)  of  several  configurations  of  an  engineering  system  for  the  removal  of  jet  fuel 
components  from  contaminated  groundwater;  these  cost  are  estimated  as  a  function  of  a  number  of 
key  engineering  and  cost  parameters.  Some  background  information  is  also  provided  to  give  the 
reader  a  brief  description  of  relevant  technical  information.  The  possible  system  configurations 
include  an  air  stripper  with  and  without  emissions  control.  Two  types  of  air  strippers  may  be  used; 
(1)  a  traditional  packed  tower  for  countercurrent  gas-liquid  contact;  and  (2)  a  centrifugal  contactor, 
also  for  countercurrent  gas-liquid  contact.  The  choice  of  gaseous  emissions  control  devices  includes: 
(1)  none,  (2)  catalytic  oxidation,  and  (3)  activated  carbon  adsorption.  Complete  computer  codes  and 
operating  instructions  for  performing  these  tasks  are  included.  The  output  from  this  assessment  is 
an  estimate  of  the  unit  processing  cost  for  remediation  of  contaminated  groundwater  for  the  selected 
engineering  and  economic  parameters. 
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SECTION  n 


RELEVANT  TECHNICAL  INFORMATION 


A.  AIRSTRIPPING 

1.  Packed  Air  Strippers 

Packed  gas-liquid  contactors  with  countercurrent  gas  and  liquid  Oows  provide  a  highly 
efficient  means  of  stripping  VOCs  from  groundwater.  The  design  of  packed  towers  is  well  developed 
in  the  chemical  engineering  literature  (References  6, 7,  and  8).  A  typical  packed  tower  is  shown  in 
Figure  3,  and  other  design  details  appear  in  Figure  4.  The  design  process  for  strippers  usually  begins 
with  known  liquid  flow  rate  and  composition  information,  including  that  of  the  solute  to  be  removed. 
The  first  step  consists  of  the  selection  of  the  stripping  medium;  in  this  section,  air  is  generally 
assumed  to  be  the  medium.  The  air  flow  rate  is  selected  such  that  an  adequate  "driving  force"  for 
this  operation  can  be  maintained.  For  conditions  that  typically  exist  in  air  strippers,  this  may  be 
satisfied  by  choosing  a  stripping  factor  greater  than  1.  The  stripping  factor  is  defined  as 

S  =  mGIL  ,  (1) 

where  m  is  the  equilibrium-phase  distribution  ratio,  for  the  VOC  of  concern.  The  stripping  factor 
also  contains  the  ratio  of  gas-to-liquid  superficial  molar  velocities,  G/L;  it  is  helpful  to  note  that  this 
ratio  is  identical  to  that  of  the  molar  flow  rates  of  these  two  phases,  since  the  tower  cross-sectional 
area  is  not  often  known  in  the  early  stages  of  the  design  process.  The  optimum  value  of  S  is 
frequently  found  to  be  between  1.25  and  4. 

Two  general  types  of  packing  are  useful  for  VOC  removal  by  air  stripping  -  random  and 
structured  -  as  shown  in  Figure  5.  Random  packings  are  available  in  a  number  of  varieties,  with 
standard  saddles  and  slotted  rings  most  commonly  used  for  commercial  applications.  The  commercial 
names  may  vary  with  the  manufacturer,  such  as  Flexirings*  from  Koch  Engineering  Company  vs  Pall 
rings  from  Chemical  Processing  Products  Division  of  Norton  Company;  and  Flexisaddles*  from  Koch 
Engineering  Company  vs  Intalox*  saddles  from  Norton  Chemical  Process  Products  or  Novalox* 
saddles  from  Jaeger  Products,  Inc.  Random  packing  is  sometimes  referred  to  as  "dumped"  packing 
due  to  the  usual  method  of  placement  in  the  tower.  Random  packings  are  available  in  a  number  of 
nominal  sizes  or  diameters,  dp,  of  up  to  0.29  feet  [89  mm]  and  in  materials  of  ceramic,  plastic,  or 
metal.  There  is  a  tendency  for  maldistribution  of  the  liquid  phase  in  towers  filled  with  dumped 
packings,  especially  when  the  packing  size  and  depth  of  packing  in  the  bed  are  incorrectly  chosen. 
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Packed  column  mteraais.  (a)  Packing  support  and  redistributor  assembly; 

(b)  Trough-type  liquid  distributor;  (c)  Perforated  pipe  distributor;  (d)  Rosette 
redistributor  for  small  towers;  (e)  Hold-down  plate,  particularly  for  low-density 
packing. 

[Reproduced  with  permission  from  Stanley  M.  Walas,  Chemical  Process 
Equipment.  (Stoneham,  MA:  Butterworth-Heinemann,  1988)] 


Figure  5.  Some  kinds  of  tower  packings:  (a)  Raschig  ring;  (b)  partition  or  Lessing  ring; 

(c)  double  spiral  ring;  (d)  metal  Flexirtng  {KOCH  Engr.  Co.)\  (e)  plastic  pall  ring; 
(f)  ceramic  Berl  saddle  {Maurice  A.  Knight  Co.);  (g)  ceramic  intalox  saddle 
{Norton  Co.);  (h)  plastic  intalox  saddle  {Norton  Co.);  (i)  metal  intalox  saddle 
{Norton  Co.);  0)  Tellerette  {Chem-Pro  Co.);  (k)  plastic  tripak  {Polymer  Piping 
and  Metals  Co.);  (1)  metal  tripak  {Pofynur  Piping  and  Metals  Co.);  (m)  wood 
grid;  (n)  section  through  expanded  metal  packing;  (o)  sections  of  expanded 
metal  packings  placed  alternatively  at  right  angles  {Denholme  Co.);  (p)  GEM 
structured  pack^g  {Glitsch  Inc,). 

[Reproduced  with  permission  from  Stanley  M.  Walas,  Chemical  Process 
Equipment  (Stoneham,  MA:  Butterworth-Heinemann,  1988}] 
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This  tendency  is  less  when  the  ratio  of  column  diameter  to  packing  size,  Dj  /dp,  is  greater  than  8 
(Reference  6),  with  the  "best"  value  occurring  at  15.  Eckert  (Reference  9)  recommends  minimum 
Di  /dp  values  of  30  for  Raschig  rings,  15  for  standard  saddles  (other  than  Raschig  rings),  and  10  to 
15  for  slotted  rings.  Because  of  the  tendency  of  the  liquid  and  gas  flows  to  segregate  in  towers  filled 
with  dumped  packings,  it  is  customary  tc  redistribute  the  liquid  at  intervals,  Zp ,  which  vary  from  2.5 
to  10  times  the  tower  diameter;  Eckert  (Reference  9)  recommends  maximum  Zp  /Dj  values  of  2.5  to 
3  for  Raschig  rings,  5  to  8  for  standard  saddles  (other  than  Raschig  rings),  5  to  10  for  slotted  rings, 
or  a  maximum  Zp  of  20  feet,  whichever  is  smaller.  There  is  still  speculation  over  the  extent  that  good 
initial  liquid  distribution  can  affect  the  maximum  Zp  /Dj  ratio  (Reference  10). 

Structured  packings  offer  advantages  of  low  pressure-drop  service  and  usually  have  excellent 
liquid  distribution  characteristics;  it  is  critical,  however,  that  the  initial  distribution  of  gases  and  liquids 
be  done  correctly  in  order  to  take  advantage  of  the  characteristics  of  these  packings.  The  structured 
packings  are  manufactured  as  elements  that  are  carefully  fitted  to  the  inside  dimensions  of  the  tower 
in  an  ordered,  or  structured,  manner.  An  example  of  such  a  packing  fabricated  from  corrugated 
sheets  is  shown  in  Figure  5.  These  sheets  are  commonly  made  of  gauze  (woven  cloth),  ceramic,  or 
sheet  metal,  as  well  as  of  various  plastics.  The  gauze  packings  are  sometimes  referred  to  as  "high- 
efficiency*  packing;  an  example  of  such  material  is  the  Koch-Sulzer  packing.  The  cloth  nature  of  the 
surface  promotes  a  capillary  action.so  that  the  liquid  covers  all  the  available  surface  even  at  low  liquid 
loadings;  the  type  of  liquid  flow  inherent  in  such  material  appears  to  greatly  enhance  liquid-phase 
transport  (References  11,  12,  and  13).  Similar  packing  elements  fabricated  of  sheet  metal  or  plastic 
do  not  appear  to  have  the  reduced  liquid-phase  resistance  and  near  constant  interfacial  area 
properties  of  the  gauze-type  packings.  The  gauze-type  structured  packing  is  more  expensive  than 
either  sheet-metal  or  plastic  structured  packing  or  dumped  packing;  its  use,  however,  can  significantly 
reduce  design  height  requirements. 

Discussion  of  typical  tower  intervals,  distributors,  packing  support,  etc.,  is  presented  by 
Treybal  (Reference  6)  and  Perry  et  al.  (Reference  14).  The  design  of  these  items  is  critical  for 
efficient  packed-tower  operation  (References  15  and  16).  They  are  usually  manufactured  and 
distributed  by  the  same  commercial  concerns  that  supply  tower  packing.  The  flow-through  types  of 
packing,  such  as  Pall*  rings,  metal  Intalox*  saddles,  and  structured  packings,  require  that  more 
attention  be  paid  to  the  distribution  of  both  gases  and  liquids  than  for  older  types  of  packings  (bluff- 
body  packings)  such  as  Raschig  rings,  Berl  saddles,  and  ceramic  Intalox  saddles.  These  older  types 
of  packing!  force  fluids  to  flow  around  them,  causing  a  higher  pressure  loss,  and  provide  more 


capacity  to  correct  for  maldistribution  (Reference  17).  Entrainment  eliminators  are  not  essential  for 
VOC  removal  in  the  air  stripper,  but  they  may  be  very  important  for  equipment  operating 
downstream  from  the  stripper  tower.  Knitted  wire  mesh  is  especially  effective  for  removing  entrained 
droplets  of  liquid  from  gas  streams,  although  many  other  devices  are  available;  further  discussion  on 
this  subject  is  found  in  References  6  and  8. 

Packed  towers  are  usually  circular  in  cross-section  due  to  ease  of  construction  and  strength 
(Reference  6).  The  diameter  of  a  tower  for  fixed  gas  and  liquid  rates  is  normally  bounded  by  limits 
of  operability.  At  a  sufficiently  small  diameter,  the  tower  will  flood.  At  too  large  a  diameter,  the 
packing  will  not  be  sufficiently  wetted  for  eflicient  mass  transfer.  Towers  have  operated  with 
superficial  liquid  velocities  as  low  as  0.007  inches/second  [0.18  mm/second];  however,  special  liquid 
distribution  systems  are  required  (Reference  18).  Usually,  the  tower  is  designed  to  operate  at  a  given 
pressure  drop  per  unit  depth  of  packing;  for  strippers,  values  of  0.25  to  0.5  inches  HjO/foot  [200  to 
400  Pascals/meter]  of  packed  depth  are  common  (Reference  6).  Alternately,  the  tower  may  be 
designed  by  selecting  a  gas  velocity  as  a  fraction  of  the  flooding  gas  velocity;  design  values  of  50  to 
80  percent  are  common. 

The  height  of  packing  required  for  a  given  separation  may  be  conveniently  calculated  using 
the  transfer  unit  concept, 

Z  =  »  ^2) 

where  Z  is  the  height  of  packing  required,  and  N^l  are  the  height  and  number  of  overall 

liquid  transfer  units,  respectively.  The  above  relationship  may  also  be  written  in  terms  of  overall  gas 
transfer  units;  the  form  chosen  usually  indicates  where  the  principal  resistance  lies.  The  height  of 
an  overall  liquid  transfer  unit  is  calculated  in  the  procedures  used  herein  from  estimates  of  heights 
of  individual  gas  and  liquid  transfer  units, 

^  Ho  /S  .  (3) 

The  correlations  for  the  individual  transfer  unit  heights  and  Hq  are  based  on  a  revised  version 
of  the  Cornell  equations,  as  presented  by  Fair  (Reference  8). 

The  equilibrium  distribution  ratio,  m,  is  related  to  Henry’s  constant,  H,  by 

m  =  y/x  =  H(p^lPt)  • 
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where  pl  ^nd  Ml  are  the  density  and  average  molecular  weight,  respectively,  of  the  liquid  phase  and 
Pt  is  the  total  pressure.. 

The  Henry’s  constant  determines  the  distribution  of  the  solute  between  the  gas  and  liquid 
phases  and  depends  upon  the  nature  of  the  solute  and  the  temperature.  Although  the  Henry’s 
constant  can  be  estimated  from  the  vapor  pressure  of  pure  solute  and  the  solubility  of  the  solute  in 
water,  experimental  determination  is  usually  recommended.  For  benzene,  the  experimentally 
determined  expression  for  the  variation  of  Henry’s  law  constant  with  temperature  that  was  used  in 
this  procedure  was  derived  from  Ashworth  et  al.  (Reference  2)  and  is  given  as 

H  =  exp(15.216  -  ft^ataj/pound  mole  ,  (5) 

where  T  is  the  absolute  temperature  in  degrees  Rankine.  Further  information  on  estimates  of  Hl 
and  Hq  can  be  found  in  a  review  by  Singh  and  Counce  (Reference  5). 

The  number  of  transfer  units  at  conditions  common  to  VOC  removal  is 


N, 


toL 


dx 

-  y/m  ’ 


(6) 


where  Xj  and  X2  are  the  liquid-phase  mole  fractions  exiting  and  entering  the  stripper,  respectively. 
The  N,ol  is  estimated  in  the  procedure  used  here  by 


to[((J^-yi/”)(i-V^/(jCt-yi/w))  >  1/5] 
(1  -  MS) 


(7) 


where  is  the  entering  gas-phase  mole  fraction  (usually  equal  to  zero).  In  its  simplest  form,  the 
number  of  transfer  units  is  the  change  in  liquid-phase  solute  composition  divided  by  the  average 
solute  driving  force  composition. 

The  design  of  packed  towers  usually  includes  a  substantial  safety  factor  to  account  for 
uncertainties  in  the  data  base.  Bolles  and  Fair  (Reference  19)  state  that  the  calculated  height  should 
be  multiplied  by  1.7  to  achieve  95  percent  confidence  when  using  their  correlations  for  H,ol. 
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The  data  base  for  estimation  of  mass  transfer  information  is  often  not  available  for  new 
packings.  Advertisement  for  new  packings  often  gives  overall  coefficient  information;  these  data, 
however,  usually  involve  the  absorption  of  CO2  into  caustic  solutions  in  which  liquid-phase  resistance 
to  mass  transfer  is  minimal  and,  thus,  is  inappropriate  for  general  use  in  the  design  of  air  strippers 
for  VOC  removal. 

In  general,  inexpensive  saddles  and  slotted  rings,  which  have  long  been  reliably  used  in 
industry,  appear  to  be  cost-effective  for  VOC  strippers  (Reference  20).  Specialty  random  packings, 
however,  are  continually  being  developed  and  marketed;  these  packings  are  usually  much  more  costly 
than  "standard"  types  of  packing.  The  economics  of  the  use  of  specialty  random  packings  should  be 
carefully  investigated  before  they  are  specified  and  used.  The  random  packing  used  for  the 
performance  and  cost  estimation  in  this  manual  is  plastic  Rexirings*.  The  major  use  of  structured 
packings  for  VOC  stripping  appears  to  be  to  repack  existing  towers  where  greater  efficiency  or  more 
capacity  is  required  (Reference  20).  Another  use  of  structured  packings  may  be  to  improve  situations 
where  minimum  space  exists  for  the  stripper  or  where  semivolatile  components  are  being  stripped. 

2.  Centrifugal  Air  Strippers 

An  alternative  to  the  traditional  packed  column  for  countercurrent  gas-liquid  contact  is  the 
centrifugal  contactor.  A  schematic  of  the  centrifugal  contactor  is  shown  in  Figure  6.  The  centrifugal 
vapor-liquid  contactor  is  composed  of  two  major  components:  the  rotating  packing  and  the  stationary 
housing.  The  liquid  phase  is  fed  into  the  center  of  the  rotating  packing  and  flows  outward  due  to 
the  centrifugal  force.  After  exiting  the  packing,  the  liquid  phase  impacts  the  housing  wall  and  flows 
by  gravity  out  of  the  unit.  The  vapor  phase  is  introduced  into  the  annular  space  between  the  packing 
and  the  housing  and  flows  inward  due  to  the  pressure  driving  force.  Seals  are  provided  between  the 
rotating  packing  and  the  housing  to  prevent  the  vapor  phase  from  bypassing  the  packing.  The  high 
shear  forces  experienced  by  the  liquid  phase  cause  the  formation  of  very  thin  films  and  rapid  renewal 
of  the  interfacial  surfaces.  The  rotation  of  the  packing  also  causes  considerable  turbulence  in  the 
vapor  phase.  Both  of  these  factors  contribute  to  efficient  mass  transfer. 

The  concepts  used  to  design  conventional  packed  towers  can  be  modified  for  the  design  of 
the  centrifugal  vapor-liquid  contactors.  This  design  is  based  on  tests  (Reference  1)  of  the  HIGEE 
centrifugal  contactor  that  is  marketed  by  Glitch,  Inc.,  of  Dallas,  Texas.  Characteristics  of  packing 
material  used  in  these  tests  are  shown  in  Table  1.  In  designing  a  conventional  packed  tower,  the 
diameter  of  the  tower  and  the  depth  of  packing  are  the  two  variables  which  need  to  be  determined. 
Similarly,  for  the  centrifugal  vapor-liquid  contactor,  the  cross-sectional  area  at  the  inner  radius  and 
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Figure  6.  A  Schematic  of  Centrifugal  Vapor-Liquid  Contactor. 
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TABLE  1.  CHARACTERISTICS  OF  PACKING  MATERIAL  USED  IN  TESTS  OF  THE 
HIGEE  CENTRIFUGAL  CONTACTOR  (REFERENCE  1) 


ROTOR 

SPECIFIC  SURFACE  AREA 
fP/ft’ 

VOIDAGE  (%) 

1 

763 

95.0 

2 

630 

93.4 

Material:  85%  nickel,  15%  chromium  (Sumitomo  Electric  Industries,  Ltd.) 

the  value  of  the  outer  radius  are  the  two  critical  variables.  An  additional  complexity  arises  in  the 
design  of  the  centrifugal  vapor-liquid  contactor  because  the  cross-sectional  area  at  the  inner  radius 
can  be  varied  by  changing  either  the  radius  or  the  axial  length.  This  results  in  an  iterative  design 
process  in  which  the  inner  radius,  outer  radius,  and  axial  length  are  varied  to  arrive  at  an  optimum 
design  solution.  The  maximum  inner  radius  is  thought  to  be  approximately  0.66  feet  [200  mm];  the 
minimum  is  0.42  feet  [127  mmj.  The  ratio  of  outer  to  inner  radius  should  be  approximately  2*. 

The  cross-sectional  area  required  at  the  inner  radius  is  dependent  upon  the  desired  hydraulic 
capacity.  This  design  of  these  units  should  be  closely  coordinated  with  the  manufacturer.  The  results 
of  hydraulic  capacity  tests  (Reference  1)  are  shown  in  Figure  7,  along  with  the  Sherwood  correlation. 
These  results  indi  te  that  the  Sherwood  correlation  underestimated  the  limit  of  operability  for  the 
rotational  speed  in  the  case  of  the  Sumitomo  packing.  However,  there  was  good  agreement  for  the 
wire  gauze  packing.  A  second-order  polynomial  curve  fit  for  the  experimental  data  is  also  shown  in 
Figure  7.  The  equation  of  this  curve  is: 

log  Y  =  -2.274484  -  1.13671og(X)  -  0.1681 18[log(X)]*  ,  (8) 

with  a  coefficient  of  determination  (r  ^ )  of  0.80.  The  Y  and  X  in  this  equation  are  the  ordinate  and 
abscissa,  respectively,  for  Figure  7. 
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*S.  P.  Singh,  Personal  Communication  with  R.  M.  Counce,  September  1989. 
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Figure  7.  Comparison  of  Limit  of  Operability  Data  with  that  Predicted  by  the  Sherwood  Flooding  Correlation. 


The  pressure  drop  across  the  rotating  packed  torus  of  a  centrifugal  vapor-liquid  contactor  is 
modeled  as  a  function  of  two  terms.  The  first  term  accounts  for  the  pressure  drop  due  to  rotation 
of  the  packing,  and  the  second  term  accounts  for  pressure  drop  resulting  from  the  flow  of  fluids 
through  a  porous  media: 

APt  =  1.356  X  10-*Pq  -  rf)  +  1.458  x  IQ-’-^Pofr^  -  atm  ,  (9) 

e 

where  pq  is  the  gas  density  in  pound/ft^  <a  is  the  rotational  speed  in  rad/second,  T2  and  rj  are  the 
outer  and  inner  radii  in  feet,  respectively,  ap  is  the  packing  specific  area  in  1/feet,  e  is  the  packing 
void  fraction,  and  is  the  average  superficial  gas  velocity  in  feet/second.  The  coefficient  of 
determination  (r  ^ )  for  the  regression  fit  is  0.94.  Although  the  approach  outlined  above  is  a  rather 
simple  representation  of  a  complicated  system,  it  does  a  reasonable  job  in  describing  the  experimental 
data  and  is  convenient  to  use. 

The  mass  transfer  results  for  the  HIGEE  centrifugal  contactor  are  based  on  a  study  by 
Wilson  et  al.  (Reference  1);  a  variation  of  the  transfer  unit  concept  was  developed  by  Singh 
(Reference  21)  to  represent  the  results  of  this  study: 

-  rf)  =  .  (10) 

where  the  area  of  an  overall  liquid-phase  transfer  unit  is 

Kl  =  (LViiL  (pWpI  s’)'®"  •  (11) 

s 

where  L*  is  liquid  mass  velocity,  pi  is  liquid  viscosity,  pl  liquid  density,  and  a^  is  acceleration.  The 
dimensionless  groups  are  the  Reynolds  number  and  Grashof  number.  Although  the  coefficient  of 
determination  (r  ^ )  from  regression  analysis  was  only  0.61,  the  correlation  predicts  the  area  of  a 
transfer  unit  (ATU)  within  plus  or  minus  20  percent,  which  is  similar  to  what  existing  correlations 
predict  for  the  conventional  packed  tower.  The  above  correlation  is  based  on  the  assumption  that 
the  rate  of  mass  transfer  is  liquid-film  controlled.  The  proposed  correlation  could  be  made  more 
general  by  including  the  Schmidt  number.  Howe\^r,  since  the  Schmidt  number  in  the  experimental 
data  remained  constant,  it  may  be  considered  to  be  part  of  the  coefficient. 

Power  consumption  for  a  centrifugal  vapor-liquid  contactor  can  be  modeled  using  two  distinct 
terms.  The  first  term  can  be  used  to  account  for  all  the  frictional  losses  and  the  second  term  to 
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account  for  the  power  required  to  accelerate  the  liquid  entering  the  packing  torus  to  the  rotational 
speed  at  the  outer  radius.  The  frictional  losses  are  highly  dependent  upon  the  design  of  the  machine 
and  cannot  be  predicted  without  advanced  knowledge  of  the  design.  The  power  required  to 
accelerate  the  liquid,  on  the  other  hand,  can  be  described  by  a  theoretical  model  (Reference  22). 
An  expression  for  the  overall  power  consumption  was  obtained  by  Wilson  et  al.  (Reference  1).  For 
the  three  rotors  tested,  the  experimental  data  for  the  region  where  the  rotational  speed  was  greater 
than  its  limit  of  operability  gave  the  following  equation: 

P'  =  1.22  +  4.64  X  10-*  Ql  . 

where  Ql  is  the  volumetric  liquid  rate.  The  coefficient  of  determination  (r  ^ )  for  this  equation  was 
0.92.  The  correlation  provides  a  reasonable  prediction  of  the  power  consumption  over  the  operating 
conditions  studied  by  Wilson  et  al.  (Reference  1). 

B.  EMISSIONS  CONTROL 

Although  air  stripping  is  an  effective  technology  for  removing  VOCs  from  groundwater,  it  simply 
transfers  contaminant  from  one  medium  to  another.  In  many  cases,  the  removal  of  the  VOCs  from 
the  air  may  be  necessary  to  reduce  the  exposure  to  these  chemicals.  Emissions  control  technologies 
have  been  used  by  industries  which  use  solvents  or  other  VOCs  in  production  processes  to  either 
recover  or  destroy  these  compounds.  The  three  most  commonly  used  techniques  for  controlling 
VOCs  from  air  streams  are  (1)  adsorption  onto  activated  carbon,  (2)  catalytic  oxidation,  and 
(3)  thermal  oxidation.  This  work  focuses  on  VOC  control  by  means  of  adsorption  onto  activated 
carbon  and  catalytic  oxidation. 

1.  Activated  Carbon 

Activated  carbon  can  be  produced  from  petroleum  fractions,  wood,  coconut  shells,  and  coal, 
and  then  given  treatment  by  superheated  steam  to  extend  the  pore  network  of  the  particles  and  to 
give  activated  carbon  a  large  surface  area. 

The  use  of  activated  carbon  to  remove  solvents  from  airstreams  is  well  established  and  has 
been  employed  since  the  1930s  (Reference  23).  The  adsorption  of  the  solvent  molecules  onto 
activated  carbon  is  mainly  due  to  Van  der  Waal  forces,  and  no  chemical  reaction  takes  place.  The 
adsorption  of  the  VOCs  from  the  gas  stream  onto  the  activated  carbon  bed  depends  upon  (1)  type 
of  carbi)n,  (2)  relative  humidity,  (3)  temperature,  (4)  concentration  and  type  of  VOC,  and 
(5)  regeneration  steps  used  (Reference  24). 
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Several  equations  are  available  in  the  literature  to  describe  the  adsorption  of  compounds  onto 
activated  carbon.  These  include  the  Braunauer-Emmett-Teller  (BET),  Freundlich,  Langmuir,  and 
Dubinin-Radushkevich  equations.  A  commonly  used  isotherm  describing  the  adsorption  of  single¬ 
component  VOCs  from  gas  streams  is  the  Freundlich  isotherm. 


(13) 


where  qj  is  the  component  loading  on  the  carbon  for  component  i,  is  the  Freundlich  coefficient, 
C;  is  the  gas-phase  concentration,  and  l/nj  is  the  Freundlich  exponent. 

Equilibrium  information  for  gaseous  VOCs  over  activated  carbon  is  difficult  to  obtain  from 
the  open  literature.  A  method  for  estimation  of  such  information  is  presented  by  Crittenden  et  al. 
(Reference  25);  this  method  is  based  on  work  by  Dubinin  and  Radushkevich  (Reference  26): 


<ii  = 


oKa 


10 


exp 


RT  In— 


(14) 


where  W<,  is  the  maximum  available  space  on  the  carbon,  is  the  liquid  density  of  the  pure 
adsorbate,  T  is  temperature  in  Kelvin,  B  is  the  microporosity  constant,  $  is  the  affinity  coefficient, 
and  Pj  and  P,  are  the  partial  and  saturation  pressures  of  component  i,  respectively.  The  parameters 
Wg  and  B  are  characteristic  of  the  adsorbent.  These  constants  for  the  reference  compound,  toluene, 
adsorbing  on  a  typical  carbon  are  imbedded  in  the  program  SPEQ,  which  is  used  to  estimate 
Freundlich  parameters.  This  program  is  included  in  this  manual  in  Appendbc  A.  The  affinity 
coefficient  can  be  determined  from  the  polarizabilities  of  the  adsorbate  and  a  reference  compound 
by: 

P  =  H/P*  . 

The  polarizability  of  the  reference  compound  is  also  imbedded  in  the  program  SPEQ.  The 
polarizability  of  any  compound  can  be  calculated  from  the  refractive  index  using  the  following 
equation: 

p  =  [(n^  -  1)M1  /  [(n2  +  2)p  J  .  (16) 

Equation  (14)  may  be  used  to  generate  qj  vs  C,  information,  which  may  be  correlated  into  the  form 
of  the  Freundlich  equation  (Reference  25). 
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Multicomponent  gaseous  VOC  adsorption  equilibria  have  been  estimated  by  Crittenden  et 
al.  (Reference  27),  using  simplified  ideal  adsorbed  solution  theory  (lAST).  This  approach  utilizes 
single-component  Freundlich  isotherm  information  to  predict  loading.  The  velocity  of  the  loading 
wave  front  for  single-component  adsorption  of  component  i  onto  activated  carbon  is  expressed  as 

V  =  -  ,  (17) 

^  (qiPc  *  Pf  Coi®) 


where  is  the  wave  front  velocity,  Vp  is  the  fluid  superficial  velocity,  is  the  carbon  bulk  density, 
C^,i  is  the  influent  concentration  of  i,  and  c  is  the  carbon  void  fraction.  This  reduces  to  the  following 
equation  for  most  cases: 


(18) 


The  wave  velocity  of  the  kth  component,  for  k  ^  2,  of  a  multicomponent  mixture  adsorbing 
onto  activated  carbon  is  given  by  Speth  (Reference  28): 


-  E  Pc  ♦  Pr  -  ViJ  Pc  *  Pr 

_ tl _ 1 _ 

(<lk.k  Pc  ^  Pf  «) 


The  subscript  i,j  indicates  compound  i  in  zone  j,  and  is  the  influent  concentration  of  component 
k;  theoretically,  there  is  a  zone  for  each  compound  in  the  multicomponent  adsorption  model.  For 
adsorption  of  N  components,  in  which  the  Nth  component  is  the  least  strongly  adsorbed. 
Equation  (19)  is  the  leading  wave  front  for  all  the  adsorbing  components  when  k  =  N.  The 
concentration  of  the  ith  species  in  zone  k  may  be  estimated  by  combining  Equation  (19)  and  lAST 
to  give: 


■  N  10, 

UiK, 


i  =  1  to  N  , 


(20) 


where  N  is  the  total  number  of  components.  Based  on  lAST,  the  ratio  of  carbon  requirements  for 
a  given  cycle  time  of  a  mixture  which  is  composed  of  components  in  addition  to  component  i  to  that 
of  a  single  component  of  concentration  Q  may  be  calculated  by  the  ratio  <f>,  where 
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(21) 


♦  =  (V,k  /  . 


Thus,  activated  carbon  requirements  for  a  multicomponent  adsorption  may  be  estimated  based  on 
single-component  isotherm  information  and  the  ratio  Component  i  should  be  one  of  the  least 
strongly  adsorbed  components,  and  its  inlet  concentration  Cj  should  be  the  same  for  both  wave 
velocity  calculations. 

The  adsorption  of  VOCs  from  the  vapor  phase  can  be  affected  by  the  relative  humidity.  This 
effect  can  be  explained  by  examining  the  adsorption  mechanism  of  water  from  the  gas  stream.  At 
low  relative  humidity,  the  removal  of  water  from  the  vapor  phase  is  due  mainly  to  adsorption  onto 
the  surface  of  the  adsorbent  Due  to  a  relatively  low  number  of  hydrophilic  sites  available  on  the 
activated  carbon  surface,  the  adsorption  of  water  is  fairly  low.  As  the  concentration  of  the  water  in 
the  vapor  phase  increases,  capillary  condensation  begins  to  occur  in  the  pores  of  the  activated  carbon 
and  the  sites  available  for  direct  adsorption  of  the  organic  molecules  from  the  vapor  phase  are 
reduced.  The  capillary  condensation  of  water  increases  the  resistance  to  mass  transfer  of  other 
components  in  the  pores  because  the  compound  must  flrst  dissolve  in  the  water  and  then  travel 
through  the  water  to  the  adsorption  site.  The  reduction  in  the  adsorption  of  VOCs  from  the  vapor 
phase  decreases  rapidly  as  the  relative  humidity  rises  above  45  to  50  percent  (Reference  29).  Thus, 
the  lifetime  of  the  adsorption  bed  can  be  extended  by  lowering  the  relative  humidity  of  the  gas 
stream,  which  can  be  accomplished  by  raising  the  gas  stream  temperature.  Increasing  the  temperature 
too  much,  however,  can  also  reduce  the  adsorption  capacity.  The  temperature  and  relative  humidity 
need  to  be  adjusted  to  obtain  the  optimum  adsorption  capacity.  For  the  calculations  in  the 
spreadsheet  described  in  this  manual,  the  gas  stream  is  heated  to  KKl^F. 

a.  Nonregieiierative 

In  many  cases,  the  amount  of  activated  carbon  used  for  vapor  phase  removal  of  VOCs  is 
small  enough  so  that  disposal  of  the  spent  carbon  is  a  viable  economic  option.  However,  other  issues 
need  to  be  considered,  such  as  potential  classification  of  the  spent  carbon  as  a  hazardous  waste.  One 
option  available  may  be  the  regeneration  of  the  carbon  by  a  commercial  vendor,  which  may  be  less 
expensive  than  on-site  regeneration. 

b.  Regenerative 

Steam  stripping  of  activated  carbon  is  a  common  method  of  reactivating  carbon  when  the 
adsorption  involves  concentrated  solvent  vapors.  Decisions  on  steam  stripping  as  a  method  of 
regeneration  are  very  dependent  on  whether  the  condensate  from  this  stripping  operation  will 
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separate  into  two  phases  for  easy  removal  of  the  organic  phase.  For  adsorption  of  dilute  VOC 
compounds  onto  activated  carbon,  the  condensate  from  steam  regeneration  may  not  contain  sufficient 
VOCs  to  form  a  separate  phase  and  may  require  treatment  of  a  secondary  aqueous  phase 
contaminated  with  dilute  VOCs.  To  compound  these  problems,  Crittenden  et  al.  (Reference  29) 
report  very  high  steam  requirements  for  regeneration  of  activated  carbon  used  in  VOC  adsorption 
service.  Some  estimated  on-site  costs  for  large-scale  regeneration  by  multihearth,  fluidized-bed,  and 
infrared  technologies  are  reported  by  Adams  et  al.  (Reference  30).  The  current  availability  of 
fluidized-bed  and  infrared  regeneration  technologies  is  questionable. 

2.  Catalytic  Qndation 

Another  type  of  technology  that  has  been  used  in  industry  to  control  the  emission  of  VOCs 
in  gas  streams  is  catalytic  oxidation.  In  this  process,  a  catalyst  is  used  to  promote  the  oxidation  of 
organic  compounds  at  lower  temperatures  than  required  for  thermal  oxidation.  The  catalyst  increases 
the  rate  of  the  reaction  by  bringing  the  reactants  together  or  by  lowering  the  activation  energy  of  the 
reaction.  Approximately  500  to  2000  catalytic  oxidation  units  are  used  to  control  the  emission  of 
VOCs  in  various  industries  (Reference  31). 

The  performance  of  a  catalytic  oxidation  unit  depends  upon  temperature,  type  and 
concentration  of  compounds,  space  velocity  (residence  time),  and  type  of  catalyst.  A  typical  catalytic 
oxidation  unit  is  usually  composed  of  four  basic  parts.  A  preheater  is  used  to  bring  the  incoming  gas 
stream  to  the  required  temperature  before  entering  the  catalyst  bed.  This  temperature  is  dependent 
on,  for  example,  the  contaminant  concentration  in  the  gas  stream,  the  type  of  catalyst  bed  (fixed-bed 
or  fluidized-bed),  and  the  operating  temperature  of  the  catalyst  bed.  After  the  preheater,  a  mixing 
chamber  is  used  to  promote  uniform  temperature  of  the  gas.  The  catalytic  system  then  follows,  which 
may  be  a  fixed-bed  or  fluidized-bed  system.  The  catalyst  is  usually  composed  of  either  metal  oxides 
or  finely  divided  precious  metal  on  either  a  metal  or  ceramic  support  structure.  The  operating 
temperature  of  the  catalyst  bed  is  determined  by  the  required  contaminant  destruction  efficiency. 
For  the  same  type  of  catalyst,  some  species  require  higher  oxidation  temperatures  than  others.  Also, 
for  the  same  contaminant  species,  the  required  oxidation  temperature  varies  with  the  type  of  catalyst 
used.  The  final  part  of  the  catalytic  oxidation  unit  is  a  heat  recovery  system,  which  is  optional. 

Spivey  et  al.  (Reference  32)  conducted  a  literature  review  on  heterogeneous  catalytic 
oxidation  of  potential  environmentally  hazardous  compounds.  They  presented  an  excellent  review 
on  the  mechanism  of  catalytic  oxidation  reactions  and  a  comparison  of  metal  oxide  and  precious-metal 
catalysts.  Listed  below  are  some  findings  presented  in  their  survey  report: 
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1.  Oxides  of  copper,  manganese,  cobalt,  chromium  and  nickel  are  the  most  active  single¬ 
metal  oxide  catalysts. 

2.  Mixed-metal  oxide  catalysts  generally  have  higher  activity  than  single-metal  oxide  catalysts. 

3.  Metal  oxide  catalysts  are  less  active  than  precious-metal  catalysts;  but  the  metal  oxide 
catalysts  are  more  resistant  to  certain  poisons,  such  as  halogens,  arsenic,  lead,  and 
phosphorus. 

Palazzolo  et  al.  (Reference  33)  studied  the  destruction  of  mixtures  of  the  halogenated 
hydrocarbons  shown  in  Table  2  using  fluidized-bcd  metal  oxide  catalysts  and  reported  that  the  overall 
destruction  efficiency  varied  from  97  to  99  percent  for  all  the  mixtures  tested.  They  also  noted  that 
the  catalyst  temperature  had  the  largest  effect  on  the  destruction  efficiency,  while  other  variables, 
such  as  mixture  composition,  air-to-fuel  ratio,  space  velocity,  and  inlet  concentration,  had  only  a 
marginal  effect  on  the  destruction  efficiency.  The  destruction  efficiei.'cy  of  tetrachloroethylene  was 
the  lowest  of  all  the  compounds  tested.  The  destruction  efficiency  across  the  preheater  ranged  from 
15  to  55  percent 

In  the  study  by  Wilson  et  al.  (Reference  1),  the  noble  metal  catalyst  was  poisoned  almost 
immediately  in  field  tests.  Sulfur  compounds  were  suspected  to  be  the  cause  of  the  poisoning  of  the 
catalyst 
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TABLE  2.  MIXTURE  COMPOSITIONS  AND  TARGET  CONCENTRATIONS  FOR 
CATALYTIC  OXIDATION  TESTS  (REFERENCE  33) 


MDCrURE 

CONCENTRATION 

MIXTURE 

DESIGNATION 

LEVEL 

COMPOUNDS 

(ppmv)* 

Mixture  1 

Baseline 

Trichloroethylene 

6.3 

1,2  dichloroethylene 

14.8 

Mixture  1 

Low 

Trichloroethylene 

1.9 

1,2  dichloroethylene 

LO 

2.9 

Mixture  2 

Baseline 

Trichloroethylene 

2.7 

Benzene 

1.5 

Ethylbenzene 

5.6 

Pentane 

11.5 

Cyclohexane 

14.1 

35.4 

Mixture  3 

Baseline 

Vinyl  chloride 

7.5 

Trichloroethylene 

1.8 

9.3 

Mixture  4 

Baseline 

1,2  dichloroethane 

10 

Trichloroethylene 

10 

1 , 1 ,2-trichloroethane 

10 

Tetrachloroethylene 

10 

40 

Mixture  4 

High 

1,2  dichloroethane 

50 

Trichloroethylene 

50 

1 , 1,2-trichloroethane 

50 

Tetrachloroethylene 

50 

200 

*ppmv  =  parts  per  million  by  volume  of  the  compound. 
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SECTION  in 


EQUIPMENT  DESIGN  AND  SIZING 

In  this  manual,  the  design  of  air  strippers  for  removal  of  soluble  jet  fuel  components  from 
groundwater  is  based  on  the  removal  of  benzene.  Benzene  was  chosen  because  it  is  an  important 
contaminant  and  one  of  the  less  volatile  of  the  soluble  jet  fuel  VOCs.  All  of  the  soluble  jet  fuel 
VOCs  of  higher  volatility  (higher  Henry’s  coefficient)  can  be  assumed  to  be  removed  in  the  stripper 
with  removal  efficiencies  at  least  that  of  benzene.  It  is  more  difficult  to  generalize  stripper  removal 
efficiencies  for  compounds  of  lesser  volatility;  one  approach  would  be  to  neglect  the  removal  of  these 
compounds  for  calculation  of  stripper  performance,  while  accounting  for  their  presence  in  any 
prediction  of  carbon  bed  loading. 

The  cost  of  traditional  packed-bed  stripper  towers  is  based  on  their  physical  size.  The 
spreadsheet  developed  for  this  manual  determines  the  size  of  the  stripper  from  specification  of 
stripping  factor,  fraction  of  gas  velocity  at  flooding,  and  component  removal  efficiency.  The  selection 
of  the  stripping  factor  determines  the  ratio  of  gas  to  liquid  flow  rates  for  this  contacting  operation. 
The  selection  of  the  appropriate  fraction  of  flooding  (ratio  of  gas  velocity  to  that  at  flooding,  both 
at  the  same  ratio  of  gas  to  liquid  flow  rates)  determines  the  tower  diameter.  The  internal 
computations  of  tower  diameters  and  pressure  drop  across  the  column  are  based  on  digitalization  of 
the  Eckert  flooding  and  pressure  drop  diagram  as  presented  by  Treybal  (Reference  6). 

The  height  of  the  overall  liquid  transfer  unit  is  based  on  calculations  of  individual  gas  and  liquid 
transfer  units.  Computation  of  the  individual  liquid  and  gas  transfer  units,  Hl  and  Hq,  is  based  on 
correlations  provided  by  Bolles  and  Fair  (Reference  19).  The  number  of  transfer  units  for  a  single 
component  is  calculated  based  on  inlet  and  outlet  mole  fractions  of  benzene. 

The  algorithm  for  estimating  the  costs  of  centrifugal  air  strippers  does  not  involve  detailed  design 
of  the  unit.  Design  of  these  units  may  be  accomplished  based  on  information  presented  earlier. 
Costing  of  these  units  in  the  developed  spreadsheet  does  not  require  detailed  design. 

For  application  to  the  cleanup  of  jet  fuel  from  groundwater,  the  design  of  fixed-bed  carbon 
adsorbers  is  based  on  that  for  single-component  adsorption  of  benzene  and  a  correction  factor  to 
account  for  decreased  single-component  bed  loading  based  on  the  presence  of  other  soluble  jet  fuel 
species.  The  correction  factor  used  in  the  spreadsheet  calculations  was  determined  from  the  carbon 
adsorption  studies  described  in  Reference  1.  The  humidity  of  the  gas  stream  is  assumed  to  have  been 
reduced  to  less  than  or  equal  to  30  percent  by  heating  prior  to  adsorption.  Values  of  Kj  and  l/n;  for 
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the  Freundlich  isotherm  are  difficult  to  establish  for  gas-phase  service,  but  may  be  estimated  using 
the  computer  program  SPEQ  developed  by  R.  D.  Cortright  under  the  direction  of  J.  C.  Crittenden 
at  Michigan  Technological  University.  The  theoretical  background  may  be  found  in  a  paper  by 
Crittenden  et  al.  (Reference  34).  A  copy  of  this  program  and  required  input  information  are 
provided  in  Appendix  A.  This  program  was  used  to  calculate  the  Freundlich  parameters  for  benzene, 
as  well  as  for  trichloroethylene  (TCE).  As  discussed  later,  the  spreadsheet  program  includes  the 
capability  for  design  for  cleanup  of  TCE-contaminated  groundwater. 

To  correct  for  the  presence  of  multiple  competing  species  (other  than  water),  a  computer 
program,  ETMOD,  is  provided  in  Appendix  B  along  with  a  brief  discussion  and  a  sample  input  file, 
ET.DAT.  This  program  was  developed  by  T.  F.  Speth,  P.  H.  Luft,  D.  W.  Hand,  and  J.  C.  Crittenden; 
the  theoretical  background  may  be  found  in  an  article  by  J.  C.  Crittenden  et  al.  (Reference  25).  The 
correction  factor  </>  (discussed  in  Section  II),  which  is  calculated  by  the  program,  accounts  for  the 
presence  of  multiple  components  and  is  used  to  estimate  the  size  of  carbon  beds  for  multicomponent 
systems.  In  the  absence  of  an  experimentally  determined  correction  factor,  such  as  used  for  benzene 
in  the  spreadsheet,  the  factor  <f>  may  be  used.  For  example,  since  the  carbon  usage  rate  is 
proportional  to  the  wave  front  velocity,  a  value  for  of  2.7029/1.3794  =  1.96  for  n-pentane  is 
calculated  from  the  output  results  in  Exhibit  B-2,  Appendix  B. 

The  algorithm  for  calculation  of  the  cost  of  catalytic  oxidation  units  is  a  function  of  gas 
throughput  alone.  The  fuel  expense  is  a  function  of  the  operating  temperature  of  the  catalytic 
oxidation  unit.  For  90  percent  destruction  of  benzene,  an  operating  temperature  of  the  catalyst  bed 
of  approximately  800®F  is  adequate  (Reference  35).  A  conservative  value  of  900 ®F  is  used  in  the 
spreadsheet  calculations.  Some  information  on  destruction  of  various  gaseous  mixtures  of  VOCs  is 
presented  by  Tichenor  and  Palazzola  (Reference  36).  Proper  selection  of  this  operating  temperature 
is  critical  to  the  achievement  of  desired  removal  efficiencies. 

The  costs  of  pumps,  as  well  as  operating  expenses,  are  related  to  the  power  requirements.  The 
total  power  requirements  for  all  pumps  supplying  groundwater  are  estimated  based  on  the  head  at 
the  pump  discharge.  This  discharge  head,  which  is  largely  due  to  the  hydrostatic  head  Zj  and  the 
head  loss  due  to  friction,  may  be  estimated  as 

=  (gig,)  Zj.  *  ,  (22) 

where  /  is  the  Fanning  friction  factor,  g  is  the  acceleration  of  gravity,  g^  is  the  gravitional  conversion 
factor,  Zt,  is  the  total  length  of  water  piping  of  diameter  D,  and  V  is  the  water  velocity  in  the  piping. 
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The  velocity  of  water  through  a  well-designed  piping  system  is  assumed  to  be  5  feet/second 
(Reference  37).  A  typical  valve  value  of  the  Fanning  friction  factor,  /,  is  0.008.  The  required  total 
power  consumption  of  the  pumps  may  be  estimated  as 

L'  z' 

p  =  -  ■  (23) 

*  3962  n, 

where  is  the  water  flow  rate  in  gallons  per  minute  and  the  efficiency  of  centrifugal  pumps  may  be 
estimated  as 

II,  =  0.7  .  (24) 

The  cost  of  the  blower,  as  well  as  operating  expenses,  is  also  related  to  the  power  consumption, 

Pg  =  G' Zg/  6350  x\g  ,  (25) 

where  G'  is  the  volumetric  air  flow  rate  in  standard  cubic  feet  per  minute,  Zg  is  the  total  blower 
static  head  in  inches  of  H2O,  and  the  overall  efficiency  may  be  estimated  (Reference  38)  as 

t1,  =  0.70.  (26) 

The  head  requirements  are  primarily  those  due  to  flow  through  the  relevant  major  equipment  items. 
The  estimation  of  the  gas-phase  pressure  loss  for  the  traditional  packed  tower  was  mentioned  earlier. 
The  pressure  drop  for  the  centrifugal  air  stripper  is  from  a  correlation  presented  earlier  in  Section  II. 
The  pressure  drops  across  the  catalytic  oxidation  unit  and  the  activated  carbon  bed  were  assumed  to 
be  5  and  4  inches  of  water,  respectively. 
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SECTION  IV 


ESTIMATION  OF  COSTS  FOR  REMOVAL  OF  VOCs  FROM  GROUNDWATER 

This  section  focuses  upon  the  estimation  of  the  cost  (in  current  year  dollars/lOOO  gallons  of  water 
processed)  of  the  removal  of  VOCs  from  groundwater.  The  estimation  procedure  involves  the 
following  sequence,  all  of  which  is  handled  on  spreadsheets: 

1.  Engineering  design  characterization,  making  use  of  design  algorithms  as  discussed  in  the 
preceding  section. 

2.  Estimation  of  (1)  the  fixed  capital  costs  for  capital  equipment  and  (2)  annual  operating  costs. 
All  costs  presented  for  this  study  are  adjusted  to  1990  dollars  using  the  estimates  for  the 
average  inflation  rate  and  the  average  annual  real  interest  rate,  as  described  later.  The  costs 
may  be  accurate  to  approximately  plus  or  minus  30  percent  overall. 

3.  Use  of  inputs  &om  these  two  steps  in  the  operating  lifetime  analysis  to  estimate  the  cleanup 
costs  per  1000  gallons  of  groundwater  processed. 

The  spreadsheet  programs  allow  for  many  of  the  engineering  and  the  financial  parameters  to  be 
treated  as  input  variables,  which,  in  turn,  allows  for  single-variable  sensitivity  analysis.  Because  of 
this,  it  is  necessary  to  establish  a  set  of  "base  case"  conditions  summarized  in  Table  3.  Some  special 
variations  in  the  base  case  parameters  are  identified  as  part  of  the  subsequent  data  analyses. 

VOC-2.100  is  a  dual  purpose  spreadsheet  that  generates  design  parameters  necessary  for 
estimation  of  fixed  capital  and  annual  operating  (noncapital)  costs.  The  design  information  is  then 
carried  over  into  estimation  (in  1990  dollars)  of  the  fixed  capital  and  the  annual  expense  items.  A 
government-owned  (break-even,  or  zero-profit)  scenario  is  stipulated.  Thus,  there  is  no  consideration 
of  the  time  phasing  between  expense  and  revenue,  and  working  capital  need  not  be  included  in  the 
cost  estimations.  The  design  for  cleanup  of  groundwater  contaminated  with  jet  fuel  components  is 
estimated  using  the  concentration  of  benzene  in  the  water  and  accounting  for  the  presence  of  other 
components  where  necessary,  as  discussed  previously.  In  addition  to  the  calculations  for  removal  of 
jet  fuel  components,  VOC-2.100  has  the  capability  to  treat  cleanup  of  TCE.  For  completeness  of 
economic  evaluation,  calculations  are  made  for  catalytic  oxidation  of  off-gas  containing  TCE. 
However,  it  is  noted  that  (1)  some  catalysts  are  poisoned  by  halogens,  and  (2)  an  environmental  issue 
exists  as  to  the  form  in  which  chlorine  is  emitted  from  an  oxidation  unit. 
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3.  INPUT  DESIGN  VARIABLES  FOR  SPREADSHEET  (TABLE  1  OF  SPREADSHEET  VOC-2.100) 


TABLE  3.  DESIGN  INPUT  VARIABLES  FOR  SPREADSHEET  (TABLE  1  OF 
SPREADSHEET  VOC-2.100)  (CONCLUDED) 


Overall  Operating  Qcle  Load  Factor 
(365  days/year  =  1.00) 

0.85 

Cost  Parameters 

Equipment  Installation  Factor  (Stripper) 

2.20 

EDF  (Carbon  and  Catalytic) 

1.60 

Materials  (SS)  Factor,  STRIPPER 

1.70 

Fuel  Oil,  1990  $/10®  BTU 

6.15 

Electricity  Cost,  1990  %/kwh 

0.060 

Operations  &  Maintenance  Factor 

0.150 

Overhead  Rate  (%)  on  Expense 

100 

Average  Annual  Inflation:  1977-90 

0.0484 

*Treybal  (1980). 


For  operating  lifetime  financial  analyses,  the  output  values  of  fixed  capital  costs  and  of  annual 
expense  from  VOC-2.100  are  used  as  inputs  to  spreadsheet  VOC-2.300.  This  operating  lifetime 
analysis  spreadsheet  provides  the  total  dollar  cost  per  1000  gallons  of  groundwater  processed  over 
a  20-year  equipment  lifetime.  The  cost  estimation  procedures  used  are  in  agreement  with  guidelines 
of  the  United  States  General  Accounting  Office  (GAO).  The  operating  lifetime  analysis  was  used 
because  the  estimation  of  capital  and  annual  operating  costs  is,  by  itself,  not  adequate  for  the 
economic  evaluation  of  a  proposed  project.  The  determination  of  total  dollar  cost  involves 
(1)  calculating  the  annual  capital  and  expense  costs  for  each  year  of  operation  and  (2)  summing  these 
yearly  costs.  The  operating  cost  for  each  year  of  the  operating  lifetime  is  calculated  using  the  annual 
expense  estimate  for  the  base  year  and  adjusting  it  to  the  current  year  with  the  assumed  inflation 
rate.  Since  the  assumed  depreciation  period  was  20  years,  the  annual  capital  cost  is  5  percent  of  the 
total  capital  costs  of  the  major  equipment  items.  The  cost  per  1000  gallons  of  groundwater  processed 
is  then  the  total  current  year  cost  (the  sum  of  the  yearly  costs)  divided  by  the  total  amount  of 
groundwater  processed  (in  thousands  of  gallons)  during  the  operating  lifetime  of  the  air  stripping 
system. 
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The  components  of  the  groundwater  cleanup  systems  consist  of: 

•  A  packed  column  stripper  or  a  centrifugal  stripper.  Included  with  this  unit  are  the  pumps 
and  piping  associated  with  wells  to  provide  access  to  the  groundwater. 

•  A  carbon  adsorption  system  for  cleanup  of  off-gas  from  the  stripper. 

•  Or,  as  an  alternative,  a  catalytic  oxidation  system  for  cleanup  of  off-gas  from  the  stripper. 

Thus,  sfac  equipment  configurations  are  possible.  Other  key  technical  design  parameters  used  for 

the  economic  studies  are  summarized  as  follows: 

•  Stainless  steel  construction  for  the  stripper  and  the  submerged  pump  for  the  groundwater 
system,  but  not  for  any  other  unit  operations. 

•  Plastic  Flexiring*  packing. 

•  Oil-fired  preheating  of  the  off-gas  entering  the  carbon  adsorption  system  to  reduce  relative 
humidity  and  improve  adsorption. 

•  Oil-Ored  heating  of  the  catalytic  oxidation  unit. 

There  are  two  important  observations  related  to  the  groundwater  cleanup  costs:  (1)  for  the 
entire  range  of  sensitivity  analyses  evaluated,  the  cleanup  costs  are  dominated  by  annual  operating 
costs  rather  than  by  Hxed  capital  costs;  and  (2)  the  cost  of  carbon  adsorption  is  sensitive  to  the 
concentration  of  contaminants  in  the  groundwater,  whereas  the  cost  of  catalytic  oxidation  is  not. 

A.  ESTD^TTON  OF  THE  nXED  CAPITAL  COST 

The  fixed  capital  costs  for  grass  roots  systems  are  estimated  by  use  of  correlation  equations  from 
the  engineering  literature,  as  summarized  in  the  following  subsections.  These  values  are  then 
escalated  by  use  of  an  equipment  installation  factor  (or  Lang  factor)  to  account  for  additional  costs 
such  as  installation,  instrumentation  and  controls,  yard  improvements,  piping,  electrical,  service 
facilities,  etc.  The  adjusted  amount  then  is  the  fixed  capital  cost  used  for  subsequent  cost  analyses. 

1.  Equipmrait  Installation  Factors 

Equipment  installation  factors  in  the  range  of  about  1  to  6  are  found  in  miscellaneous  sources 
(References  37, 39,  and  40).  Explicit  values  for  selected  pollution  control  equipment,  however,  seem 
to  lie  in  the  range  of  about  1.6  to  2.5.  The  base  case  values  used  in  this  study  were  2.0  for  either  of 
the  two  strippers,  and  1.6  for  either  of  the  two  air  treatment  systems.  Either  or  both  may,  however, 
be  used  as  input  variables. 

2.  Fned  Capital  Investment  for  Air  Stripper  and  Associated  Equipment 
a.  Stripper  Air  Fan  and  Motor 

Cost  data  for  air  fans  from  Vatavuk  &  Neveril  (Reference  35)  were  examined.  The  cost 
data,  however,  are  presented  in  the  fonn  of  somewhat  complex  combinations  of  graphs  which  are 
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difficult  to  translate  into  computer  algorithms.  Therefore,  cost  values  were  obtained  manually  from 
their  graphs  as  a  function  of  the  range  of  stripper  sizes  used  in  this  manual.  Since  the  stripper  size 
is  characterized  m  part  by  water  throughput  rates,  it  was  possible  to  generate  the  following 
correlation: 

C,  =  *  1000  *  [L  V301®^  _  (27) 

where  C,  is  the  cost  in  1977  dollars  of  the  air  fan  and  motor,  Fm  is  the  applicable  materials  factor, 
and  L'  is  the  water  flow  rate  in  gallons  per  minute. 

b.  Water  Pump  and  Motor 

Corripio  et  al.  (Reference  41)  use  multiple  equation  sets  for  cost  estimation  purposes; 
these  multiple  sets  are  somewhat  awkward  to  install  into  a  spreadsheet.  A  few  cost  values  were 
manually  calculated  from  their  sets  of  equations.  A  log-log  scaling  plot  was  then  evaluated,  and  the 
following  cost  correlation  equation  was  generated  for  a  submersible  centrifugal  pump  with  a  stainless 
steel  head: 

C,  =  1.7  •  4040  (171440)®^  ,  (28) 

where  C,  is  the  cost  of  pump  and  motor  in  1977  dollars. 

The  scaling  factor  of  0.7  is  consistent  with  median  values  for  many  equipment  systems 
(Reference  37).  In  the  absence  of  any  data  for  the  scaling  factor,  a  common  approach  is  to  estimate 
on  the  basis  of  0.6. 

c.  Stripper  Shell 

The  weight,  W„  of  the  total  stripper  shell  and  structural  accessories  was  estimated  on  the 
assumption  of  0.2S-inch  steel  construction. 

The  C,  (Reference  42)  in  1977  dollars  was  then  estimated  as: 

C,  =  .  (29) 


where 

Fm  =  1.0  for  carbon  steel,  1.7  for  stainless  steel, 

Cb  =  exp{6.329  +  0.18255(lnW,)  +  0.02297(lnW,)2], 
W.  =  1.1  irDjZ'T.p^. 

d.  Plastic  Pall  Ring^ 

The  cost  of  plastic  Pall  rings  in  1990  dollars  is: 
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Cf  =  %l\.QOIft^  {1-inch  plastic  Floorings) 


(30) 


Cf  =  %5.50lft^  (2-inch  plastic  Flexirings)  (31) 

=  $3.50//?’  (3 -inch  plastic  Flexirings)  (32) 

3.  Used  Cost  of  Emission  Control  Equipment  Utili^g  Activated  Carbon  Adsorption 

The  size  chosen  for  the  carbon  adsorption  units  affects  the  initial  capital  costs,  the  recycle 
interval  (the  length  of  time  between  regenerations),  and  the  operating  expenses.  In  order  to  size  the 
carbon  units,  the  number  of  days  in  a  recycle  interval  was  set  up  as  an  input  variable  in  the 
spreadsheet.  This  allows  calculation  of  a  carbon  unit  size  for  whatever  other  technical  specifications 
may  have  been  input  for  costing  purposes.  One  additional  carbon  unit  is  specified  for  standby  and 
to  minimize  downtime  during  recycle.  The  initial  purchase  requirement  for  the  carbon  load  is  then 
determined  by  the  size  and  number  of  carbon  adsorption  units. 

The  cost  of  the  carbon  units  for  the  basic  unit,  fans,  motors,  and  controls,  exclusive  of  the 
initial  carbon  load  itself,  (Reference  43)  is 

C,  »  F^(3.32  *  36,400)  ,  .  (33) 


where 

C,  is  the  cost  in  1977  dollars, 

W.  =  2 

10,000  ^  W.  ^  200,000, 

Fm  =  1.0  (carbon  steel). 

The  cost  of  the  initial  carbon  bed  load  is  based  on  the  1977  cost  of  carbon  at  $1.00/pound. 

4.  Imied  Costs  of  Emission  Control  Equipment  Utilizing  a  Catalytic  Oxidation  Unit 

Vatavuk  and  Neveril  (Reference  35)  provide  a  nonlinear  correlation  equation  for  catalytic 
oxidation  units.  In  the  size  range  of  interest  to  this  study,  however,  a  good  linear  fit  can  be  made  to 
that  portion  of  their  Figure  4  of  interest.  For  "packaged  units"  (which  includes  preheater,  blower, 
etc.),  the  cost  is 

C,  =  2.25  G'  ♦  24,000  ,  (34) 

where  G'  is  the  air  flow  rate  in  standard  cubic  feet  per  minute  and  C,  is  in  1979  dollars. 
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S.  I%GBd  Capital  Investment  Cost  for  Centrifugal  Contactor 

Actual  purchase  prices  for  two  centrifugal  contactor  units  were  available  for  50  and 
1000  gallon  per  minute  capacities*.  Log-log  scaling  gave  the  centrifugal  contactor  cost  as 

C,  *  189,000  «?j^/2.23)®^‘  .  (35) 

where  Ql  is  the  liquid  flow  rate  in  cubic  feet  per  second  and  C,  is  in  1990  dollars.  This  equation 
holds  for  0.11  Ql  ^  2.23.  The  cost  figures  supplied  were  for  complete  systems,  including  packing. 
B.  ESTIMATION  OF  OPERATING  EXPENSE  COSTS 

1.  General  Assumptions  for  Inflation  and  Intoest  Rates 

The  average  annual  inflation  rate  for  the  1977-1988  period  in  the  chemicals  industries  has 
been  4.84  percent.  This  is  a  compound  average  of  the  Chemical  Engineering  and  the  Marshall  & 
Swift  annual  inflation  indices  (provided  monthly  in  Chemical  Engineering  journal).  In  all  cases,  that 
annual  rate  has  been  used  to  adjust  older  price  bases  forward  to  1990  dollars. 

Annual  interest  generally  is  related  to  annual  inflation.  Twenty  or  more  years  ago,  applied 
interest  rates  tended  to  be  approximately  3  percent  higher  than  the  inflation  rate;  i.e.,  the  "real 
interest”  rate  was  averaging  around  3  percent  annually.  During  the  late  1970s,  the  real  interest  rate 
rose  substantially.  Over  the  past  decade,  it  has  been  moving  generally  back  down  to  the  vicinity  of 
5  percent  (Reference  44).  A  5  percent  annual  real  interest  rate  was  used  throughout  this  study, 
although  it  is  treated  as  an  input  variable  and  other  values  can  also  be  examined.  The  combination 
of  the  real  interest  rate  and  the  inflation  rate  leads  to  an  applied  interest  rate  of  slightly  less  than 
10  percent  aimually  as  the  base  case. 

2.  General  Assumptions  for  Overhead  Rates 

Cost  estimation  for  overhead  rates  may  well  generate  more  diverse  opinion  than  any  other. 
Both  a  single  rate/single  base  (References  44,  45,  46,  and  47)  and  a  multiple  rate/multiple  base 
method  were  evaluated.  The  maximum  spread  between  the  two  methods  was  about  plus  or  minus 
10  percent,  although  the  overall  uncertainty  of  the  total  analysis  may  well  be  plus  or  minus  25  percent 
(Reference  37).  Thus,  it  was  felt  that  the  more  complex  multiple-rate  method  was  not  warranted  and 
the  bulk  of  the  evaluation  was  done  with  the  single  rate/single  base  method. 


*S.  P.  Singh,  Personal  Communication  with  R.  M.  Counce,  September  1989. 
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In  general,  a  high  overhead  of  100  percent,  shown  as  an  input  variable  of  Table  3,  was  used 
for  either  the  packed  column  stripper  or  the  centrifugal  contactor.  This  was  done  because,  in  all 
cases,  all  of  the  direct  labor  is  loaded  into  the  expense  base  for  the  stripper  unit,  as  discussed  below. 
A  lower  rate  of  25  percent  was  applied  to  the  two  air  treatment  systems. 

3.  Fuel  and  Electridty  Costs 

The  annual  fuel  and  electricity  requirements  are  provided  by  the  design  portion  of  the 
spreadsheet. 

The  fuel  costs  involve  the  assumption  that  the  air  stream  will  be  preheated  using  No.  2  oil 
as  the  energy  source.  For  the  carbon  adsorption  unit,  the  air  is  preheated  to  100°F.  This  preheating 
is  to  reduce  the  relative  humidity  of  the  stream  and  thereby  improve  the  performance  of  the  carbon 
adsorption  system.  For  the  catalytic  oxidation  unit,  the  air  is  preheated  to  900 °F. 

Both  fuel  and  electrical  power  costs  vary  around  the  country.  The  base-case  values  used  here 
are  typical  of  current  costs  in  the  southeastern  United  States  in  1990;  however,  they  are  provided  as 
input  variables. 

4.  Charcoal  Regeneration 

The  regeneration  of  charcoal  is  a  major  component  of  the  annual  expense  base  for  the 
carbon  adsorption  units.  This  expense  may  be  estimated  by 


Eg  =  0.77  -  1.94 


[10) 


*  2.81 


(  wi 


[10 


-  1.38 


(36) 


where  Er  is  in  dollars  per  pound  of  carbon  regenerated.  This  correlation  expression  was  generated 
by  curve  Gtting  to  the  graphical  information  of  Adams  et  al.  (Reference  30)  for  the  multihearth 
method  of  carbon  regeneration.  The  authors  indicate  that  the  regeneration  cost  estimates  allow  for 
25  percent  excess  capacity,  about  12  percent  carbon  replacement  during  each  regeneration,  and  for 
ail  other  related  costs. 

5.  Maintenance 

Annual  maintenance  costs  in  the  range  5  to  15  percent  of  the  fixed  capital  investment  are 
commonly  quoted  (References  37  and  39).  This  wide  range  is  simply  indicative  of  the  actual  wide 
range  of  true  maintenance  costs. 

Because  of  such  uncertainties,  the  maintenance  cost  is  handled  as  an  input  variable  in 
Table  3.  A  typical  value  is  the  median  of  10  percent. 
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6.  Pall  Ring  "Maintenance” 

Total  replacement  on  a  4-year  cycle  was  assumed,  which  leads  to  an  annual  replacement  cost 
equal  to  25  percent  of  the  initial  investment  cost. 

7.  Labor 

One  full-time  equivalent  (FTE)  engineering  technician  level  at  $24,000/^ear  direct  salary  in 
1990  dollars  was  assumed.  This  is  somewhat  higher  than  suggested  by  Neveril  (Reference  39),  but 
the  higher  work  load  was  used  here  as  a  conservative  approach.  No  labor  is  loaded  onto  the  two  air 
treatment  systems.  All  labor  is  loaded  onto  the  stripper  for  all  technical  options. 

C  TECHNICAL  COMBINATIONS  OF  EQUIPMENT 

Six  technical  configurations  were  considered  in  this  study,  as  shown  below: 


Packed-column  stripper 

ALONE 

PLUS  carbon  adsorption  treatment 

PLUS  catalytic  oxidation 

Centrifugal  Contactor 

ALONE 

PLUS  carbon  adsorption  treatment 

PLUS  catalytic  oxidation 

Both  fixed  capital  and  expense  components  of  the  processing  costs  were  presumed  to  be  linearly 
additive  for  any  combination  of  the  applicable  unit  operations.  More  information  on  capital  and 
operating  costs  for  these  options  is  found  in  Wilson  et  al.  (Reference  1). 

D.  ESTIMATION  OF  CLEANUP  COSTS  BY  OPERATING  LIFETIME  FINANCIAL 
ANALYSIS 

As  described  above,  the  capital  and  annual  operating  expense  costs  from  VOC-2.100  are  input 
into  VOC-2.300  to  calculate  the  total  cost  per  1000  gallons  of  groundwater  processed  over  a  20-year 
depreciation  period.  The  operating  period  is  also  assumed  to  be  20  years  in  the  example  results 
presented  later.  By  changing  inputs  to  the  spreadsheet  (per  the  instructions  in  Section  V),  the  effects 
of  independently  varying  the  depreciation  period  and  the  operating  period  may  be  examined.  For 


I 
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example,  the  operating  period  could  be  5  years  with  the  same  depreciation  period  of  20  years.  The 
processing  costs  will  be  seen  to  decrease,  primarily  because  the  effect  of  inflation  on  expense  costs 
is  less  for  the  shorter  time.  If  the  depreciation  and  the  operating  periods  are  set  to  the  same  value, 
the  capital  cost  will  become  a  greater  fraction  of  the  total  cost  as  the  depreciation/operating  period 
decreases. 
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SECTION  V 


INSTRUCnONS  FOR  USE  OF  COST-ESTIMATING  SOFTWARE 

This  section  presents  the  instructions  for  using  the  cost-estimating  software  (available  from 
Tyndall  AFB,  as  detailed  later).  One  software  program  may  be  used  to  generate  capital  and 
operating  expense  cost  estimates  for  air  stripping  with  emissions  control  processes.  These  estimates 
may  then  be  input  into  the  operating  lifetime  financial  analysis  program  to  produce  estimates  of  the 
cleanup  costs  per  1000  gallons  of  groundwater  processed. 

In  the  following  operational/procedural  instructions,  it  is  assumed  that  the  user  is  generally 
familiar  with  spreadsheets,  and  specifically  familiar  with  the  MS  Multiplan  series  of  spreadsheets.  For 
more  explicit  information  on  the  spreadsheet  software,  refer  to  user  manuals  or  to  third  party  "user 
texts." 

PRECAUTIONARY  NOTE:  The  spreadsheets  written  for  this  project  are  in  MS  Multiplan 
Version  2.01.  They  may  be  "upward  loaded"  into  Version  3.0  or  4.0,  simply  by  using  the  TRANSFER 
LOAD  sequence.  In  the  process  of  upward  loading,  however,  the  spreadsheet  will  be  rewritten  to 
the  disk  in  the  higher  version  number  and  the  lower  version  number  will  be  lost.  Thus,  before 
executing  an  upward  -load,  a  backup  copy  in  the  onginal  version  number  should  first  be  made.  The 
spreadsheet  may  also  be  exported  into  MS  EXCEL.  Refer  to  the  MS  EXCEL  user  manual  for 
proper  instructions.  The  spreadsheet  may  be  exported  into  LOTUS  123;  however,  significant  portions 
of  the  math  and  logic  structure  may  be  lost  in  the  process.  LOTUS  and  Multiplan  do  not  use  fully 
compatible  logic,  math,  and  naming  protocols,  and  major  portions  of  these  may  have  to  be  rewritten 
if  transferred  to  LOTUS. 

The  design  and  cost-estimating  routines,  as  presented  in  the  following  Sections  A  and  B, 
generate  estimates  of  the  fixed  capital  costs  and  of  the  annual  operating  expense  costs  in  terms  of 
the  start-up  year.  These  data  are  a  necessary  but  not  a  sufficient  set  of  information  for  system 
evaluation.  For  proper  comparison  of  alternative  systems,  it  is  necessary  to  make  use  of  operating 
lifetime  analysis  informatioii,  as  in  Section  C 

A.  DESIGN  AND  COST  ESTIMATING,  WITHOUT  SENSITIVITY  ANALYSIS 

A  printout  of  the  design  and  cost-estimating  spreadsheet  program,  VOC-2.100,  is  shown  in 
Appendix  D.  The  variables  in  the  spreadsheet  are  identified  in  Appendbc  C. 

The  first  12  tables  of  VOC-2.100  are  set  up  to  handle  requisite  design  calculations.  Seven  tables, 
D-DIB  through  D-D7B,  handle  the  cost  estimating  and  summaries  for  the  jet  fuel  spill  case.  As 
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When  RECALC  is  set  to  NO,  ITERATION  can  be  set  to  YES  or  NO. 

To  locate  any  desired  input  variables  of  Table  D-1  (or  any  named  variables,  ranges,  or  areas 
anywhere  in  the  spreadsheet),  use  the  GOTO  NAME  sequence.  Using  this  sequence,  set  the  variable 
[VNUM]  equal  to  zero  to  suppress  the  sensitivity  analysis: 


1  GOTO 

NAME 

HIRSSBIHii 

0 

Again  using  the  GOTO  NAME  sequence,  enter  the  appropriate  alpha  or  numeric  information 
for  the  variables  of  interest  in  Table  D-1.  For  example. 

1  GOTO 

NAME 

1000 

will  set  the  groundv 

of  interest,  make  c 

information  in  each 

variable  has  been 

spreadsheet  may  no 

^ater  flowrate  to  1000  gallons  per  minute.  After  setting  values  for  the  variables 
:ertain  that  no  variable  has  been  set  to  [VAR].  Do  this  by  checking  the 
cell  in  each  of  the  two  columns  with  the  heading  VALUE  in  Table  D-1.  If  a 
set  to  [VAR],  enter  the  appropriate  alpha  or  numeric  information.  The 
w  be  run  by  resetting  RECALC  and  ITERATION: 

1  OPTIONS 

RECALC 

YES 

ITERATION 

YES 
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The  spreadsheet  will  run  and  generate  a  single  set  of  outputs  for  the  set  of  input  variables  as 
specified  in  Table  D-1.  Table  D-2  summarizes  those  inputs  plus  some  additional  "derived  values." 

Any  or  all  of  the  26  tables  of  the  spreadsheet,  any  range,  or  any  single  variable  value  may  then 
be  printed  by 


PRINT 


OPTIONS 


with  the  table,  range,  or  variable  specified  by  NAME  and  with  other  parameters,  such  as  margins, 
etc.,  specified  through  the  spreadsheet  print  menu. 

B.  DESIGN  AND  CX)ST  ESTIMATING,  WITH  SENSITIVITY  ANALYSIS 

The  spreadsheet  is  set  up  to  provide  the  option  of  automatic  iteration  through  any  specified 
sequence  of  values  of  any  single  numeric  input  variable  found  in  Table  D-1.  Because  of  the  logic 
structure  involved  in  "table  building"  in  this  case,  the  run  time  will  be  approximately  20  minutes  for 
an  IBM  XT-generation  machine.  The  run  time  will  be  approximately  1.5  minutes  for  a 
386-generation  machine. 

To  perform  a  single-variable  sensitivity  analysis,  proceed  as  already  outlined  above.  However, 
before  initiating  calculations  by  resetting  RECALC  to  YES,  the  following  additional  steps  must  be 
carried  out.  First,  go  to  the  variable  [VI]  using  the  GOTO  sequence: 


GOTO 


NAME 


Enter,  as  an  alpha  entry,  using  the  naming  notation  of  the  Table  D-1,  the  NAME  of  the  input 
variable  for  which  the  sensitivity  analysis  is  to  be  performed.  For  example,  if  the  groundwater 
flowrate  is  the  variable  of  interest,  then: 


ALPHA 


[WGM] 


Enter  numeric  values  for  the  next  three  variables,  [VIN],  [VSTEP],  and  [VNUM].  These  entries 
will  then  control: 
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•  the  initial  numeric  value  of  the  variable, 

•  the  incremental  value  for  each  iteration  of  the  sheet,  and 

•  the  total  number  of  values  to  be  run  (maximum  number  for  [VNUM]  is  15). 

Finally,  before  initiating  calculations,  the  variable  for  which  the  sensitivity  analysis  is  being  done 
must  be  set  to  [VAR]  using  the  VALUE  command.  For  example,  if  [WGM]  is  the  variable  of 
interest: 


GOTO 

NAME 

[WGM] 

VALUE 

mmm 

As  in  the  previous  section,  make  certain  that  all  other  input  variables  in  Table  D-1  have  been  set  to 
an  appropriate  numeric  or  alpha  value,  and  NOT  to  [VAR].  When  all  variable  entries  are  complete, 
initiate  calculations  by: 


OPTIONS 

RECALC 

YES 

ITERATION 

YES 

The  14  cost  tables  contained  in  Appendix  D  contain  an  illustrative  sensitivity  run.  These  tables 
will  automatically  clear  and  reset  themselves  if  one  subsequently  wishes  to  run  a  different  sensitivity 
analysis. 

C  OPERATING  LIFETIME  FINANCIAL  ANALYSIS 

A  complete  illustrative  printout  of  the  program  spreadsheet  VOC-2.3(X)  is  provided  as 
Appendix  E.  This  spreadsheet  is  used  for  operating  lifetime  financial  analyses.  General  procedures 
(i.e.,  the  initial  suppression  of  recalculations,  subsequent  printing)  are  the  same  as  outlined  previously. 

For  calculation  of  the  processing  cost  per  1000  gallons  of  groundwater  with  VOC-2.300,  the 
values  for  capital  and  operating  expense  costs  obtained  from  running  VOC-2.100  are  used.  As 
described  previously,  running  VOC-1.00  without  sensitivity  analysis  will  generate  one  set  of  capital 
and  operating  expense  costs,  while  running  with  sensitivity  analysis  can  generate  up  to  15  sets  of  costs. 

To  run  VOC-2.300  for  a  particular  set  of  capital  and  operating  expense  values,  first  set  RECALC 
to  NO.  Then,  in  Table  E-1,  enter  the  value  for  the  capital  cost  by: 


GOTO 

NAME 

HHuHH 

CAP 

where  CAP  is  the  numeric  value  for  capital  cost.  Next,  enter  the  operating  expense  cost  by: 
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where  EXP  is  the  numeric  value  for  the  initial  year  operating  expense  cost.  Finally,  enter  in 
Table  E-1  any  other  numerical  changes/inputs  that  you  wish  to  have  apply  to  the  equipment  system 
construction  phase. 

In  Table  E-4,  set  the  groundwater  pumping  rate  (in  gallons  per  minute)  to  the  value  that 
corresponds  to  the  set  of  capital  and  operating  expense  costs  that  was  input  in  Table  E-1: 


GOTO 

NAME 

rWPRl 

RATE 

where  RATE  is  the  numeric  value  of  the  pumping  rate.  The  variable  [M],  the  number  of  operating 
years,  may  be  set  to  some  other  value  less  than  20  years  (the  spreadsheet  is  set  up  for  a  maximum 
of  20  years),  if  desired.  Values  for  other  variables  in  Table  E-4  may  be  entered,  except  for  those 
values  which  are  derived  from  previous  inputs.  The  spreadsheet  may  now  be  run  by: 


OPTIONS 

RECALC 

YES 

ITERATION 

NO 

since  no  iteration  is  required  for  the  calculations  of  VOC-2.300. 

When  calculations  are  completed,  the  processing  cost  information  may  be  obtained  from 
Table  E-8  of  the  spreadsheet.  The  total  processing  cost  per  1000  gallons  on  a  cashflow  basis  is  found 
in  Row  8,  Column  68.  On  a  depreciation  basis,  the  total  processing  cost  is  found  in  Row  19,  Column 
68.  The  individual  contributions  to  the  total  cost  of  the  operating  expense  and  the  fixed  capital 
depreciation  are  found  in  Row  13,  Column  68  and  Row  16,  Column  68,  respectively. 

D.  EXAMPLES  OF  DESIGN  AND  COST  ANALYSIS  RESULTS 

The  information  presented  in  this  section  was  generated  by  (1)  running  spreadsheet  VOC-2.100 
to  perform  design  and  cost  estimation  with  sensitivity  analysis  to  produce  estimates  of  initial  Fixed 
capital  cost  and  initial  year  expense  cost  and  (2)  running  spreadsheet  VOC-2.3(X),  using  the  estimated 
costs  from  Step  (1)  as  input,  to  produce  the  operating  lifetime  cost  estimates.  In  Step  (1),  sensitivity 
analyses  were  performed  for  three  variables;  groundwater  flow  rate,  stripping  factor,  and  flooding 
factor.  Air  stripping  of  TCE  (trichloroethylene)  from  groundwater  was  considered,  where  the  base- 
case  conditions  were  (1)  99  percent  removal  efficiency,  (2)  a  stripping  factor  of  4,  (3)  a  flooding 
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factor  of  0.4,  and  (4)  10  ppm  of  TCE  in  the  groundwater.  Figures  8  through  10  show  the  dollar  cost 
per  1,000  gallons  of  groundwater  processed  as  a  function  of  groundwater  flow  rate,  stripping  factor, 
and  flooding  factor,  respectively.  Included  on  each  figure  are  the  results  for  the  three  different  cases 
of  (1)  air  stripping  only,  (2)  air  stripping  with  carbon  adsorption,  and  (3)  air  stripping  with  catalytic 
oxidation. 

Figure  8  demonstrates  the  economy  of  scale;  the  processing  cost  decreases  with  increasing  water 
flow  rate  or,  in  effect,  the  size  of  the  equipment  There  is  little  difference  between  the  packed 
stripper  and  the  centrifugal  contactor  costs.  As  expected,  the  costs  for  the  two  cases  with  emissions 
control  are  significantly  higher  than  for  air  stripping  alone.  For  TCE,  the  carbon  adsorption  and  the 
catalytic  oxidation  emissions  control  options  produce  approximately  the  same  processing  costs  as  a 
function  of  water  flow  rate.  With  higher  concentrations  of  contaminant  in  the  groundwater,  the 
carbon  cost  will  increase  due  to  the  higher  carbon  usage  while  the  catalytic  oxidation  cost  will  be 
essentially  unaffected.  This  may  be  seen  in  Wilson  et  al.  (Reference  1),  where  cases  are  presented 
for  different  concentrations  of  jet  fuel  in  groundwater. 

In  Figure  9,  the  processing  cost  for  air  stripping  alone  is  relatively  insensitive  to  stripping  factor. 
However,  for  carbon  adsorption  and  catalytic  oxidation,  the  cost  increases  with  increasing  stripping 
factor.  In  the  case  of  catalytic  oxidation,  the  cause  of  the  increasing  cost  is  the  larger  size  of  the 
oxidation  unit  required  as  the  air  flow  rate  increases.  For  carbon  adsorption,  the  capacity  of  the  bed 
decreases  with  increasing  air  flow  rate  as  a  result  of  the  lower  concentration  of  the  contaminants  in 
the  effluent  air  from  the  stripper.  Figure  10  shows  that,  in  all  three  cases,  the  processing  costs  are 
independent  of  flooding  factor  for  all  practical  purposes. 

As  pointed  out  in  Section  V.A,  there  are  39  input  variables  to  the  cost  estimation  spreadsheet 
which  may  be  examined  for  their  impact  upon  the  cost  of  the  air  stripping  systems.  Examples  of  the 
effects  of  three  of  these  variables  have  been  presented  here.  In  Appendix  E,  results  from  single 
variable  sensitivity  analyses  are  given  for  installation  factor,  contaminant  removal  efficiency,  the 
applied  overhead  rate,  and  the  carbon  regeneration  interval. 

The  spreadsheets  described  here  should  be  very  useful  for  the  evaluation  of  competing 
alternatives  for  groundwater  cleanup.  Quite  often,  the  choice  between  alternatives  is  made  on  the 
basis  of  competitive  bids  for  the  installed  capital  equipment  system,  i.e.,  the  low  bidder  on  equipment. 
However,  as  may  been  seen  in  the  detailed  analyses  presented  in  the  Appendix,  one  should  not 
overlook  the  significant  effect  of  the  operating  exp>ense  in  determining  the  true  cost  of  cleanup 
operations. 
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Rgure  8.  Processing  Cost  vs  Water  Flow  Rate  for  TCE  Cleanup. 
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Copies  of  the  VOC-2.100  and  VOC-2.300  spreadsheets  are  available  from  the  Chemical/Physical 
Treatment  Technology  Area  Manager,  HQ  AFESC/RDVW,  Tyndall  AFB,  FL.  Include  a  diskette 
with  your  request. 


47 


SECTION  VI 


CONCLUSIONS 

The  spreadsheet-based  mcxlel  developed  and  described  herein  can  be  readily  utilized  by  Air 
Force  personnel  for  the  evaluation  and  selection  of  an  air  stripping  with  emissions  control  technology 
for  the  remediation  of  VOC  contaminated  groundwater.  If  desired,  the  correlations  for  design  and 
cost  estimations  can  be  updated  as  new  information  becomes  available.  By  varying  financial  input 
parameters,  such  as  depreciation  period,  interest  rates,  etc.,  different  scenarios  may  be  evaluated  for 
particular  operating  conditions  of  groundwater  flowrate,  contaminant  concentration,  etc. 
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APPENDIX  A 


ESTIMATION  OF  GAS-PHASE  FREUNDUCH  ISOTHERM  PARAMETERS 
FOR  VOCS  OVER  ACTIVATED  CARBON 

To  estimate  the  performance  of  gas-phase  adsorption,  adsorption  equilibria  are  necessary.  The 
program  SPEQ,  furnished  by  J.  C.  Crittenden  of  Michigan  Technological  University,  provides 
estimates  of  such  equilibria  in  the  form  of  single-component  Freundlich  parameters.  Application  of 
this  program  is  presented  by  Crittenden  et  al.  (Reference  27).  The  program  utilizes  the  Dubinin- 
Radushkevich  equation  described  in  Chapter  2.  The  program  proceeds  interactively  and  requires  the 
following  information  for  each  component  of  interest: 

Compound  name 
Antoine  Equation  Constant  A 
Antoine  Equation  Constant  B 
Antoine  Equation  Constant  C 
BETA  (affinity  coefficient) 

Density  (gram/cm^) 

Molecular  Weight 
System  Pressure  (mm  Hg) 

Temperature  ("C) 

Desired  Gas  Concentration  (micrograms/liter  at  STP) 

Initial  concentration  for  integration  (;igram/liter) 

Integration  concentration  step  (/tgram/liter) 

The  latter  two  parameters  must  sometimes  be  determined  by  trial  and  error,  but  may  be  set  at 
about  1/100  and  1/1000,  respectively,  of  the  desired  gas  concentration.  The  program  output  is  found 
in  SPEQ.OUT.  This  output  includes  the  Freundlich  parameters  K  [(/imol/gram)/(/zmol/L)‘^]  and  1/n, 
as  well  as  the  gas  concentration  range  for  plus  or  minus  10  percent  prediction  accuracy  loading.  The 
program  and  sample  output  are  shown  as  Exhibits  Al  and  A2,  respectively. 
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EXHIBIT  A-1 


prcx;ram  name  -  spbq 

LANGUAGE  -  FORTRAN 

THIS  PROGRAM  USES  THE  D-R  EQUATION  TO  PREDICT  THE  SURFACE  LOADING  FOR 
A  DESIRE  GAS  CONCENTRATION  (ug/L  STP) .  THE  PROGRAM  WILL  THEN  FIND 
THE  CORRECT  FREUNDLICH  PARAMETER  IS  SUCH  A  MANNERS  AS  THE  SPREADING 
PRESSURE  FOR  THE  CALCULATED  FREUNDLICH  PARAMETER  WILL  BE  THE  SAME  AS 
THE  SPREADING  PRESSURE  CALCULATED  FOR  THE  D-R  EQUATION. 

THE  SPREADING  PRESSURE  IS  FOUND  BY  DOING  A  NUMERICAL  INTEGRATION  THAT 
MAKES  USE  OF  THE  TRAPEZOID  RULE. 

PROGRAM  DEVELOPED  BY:  RANDY  D.  CORTRIGHT,  Graduate  Student 

DEVELOPED  AT  :  MICHIGAN  TECH  UNIVERSITY 

Houghton,  Michigan  49931 


INPUT: 

INPUT  THE  FOLLOWING  PARAMETER  INTO  A  FILE  NAME  SPEQ.IN 

COMPOUND  NAME 
ANTOINES  EQTN.  A 
ANTOINES  EQTN.  B 
ANTOINES  EQTN.  C 
BETA 

DENSITY  G/CC 
MOLECULAR  WEIGHT 
SYSTEM  PRESSURE  IN  MM  HG 
TEMPERATURE  DEG  C 

DESIRE  GAS  CONCENTRATION  IN  MICROGRAMS/LITER  (STP) 


OUTPUT 

OUTPUT  WILL  BE  SENT  TO  THE  FILE  NAME  SPEQ.OUT 

THE  OUTPUT  WILL  GIVE  INTERMEDIATE  SPREADING  PRESSURE  VALUES  FOR 
INTERMEDIATE  GAS  CONCENTRATIONS,  THE  FINAL  SPREADING  PRESSURE, 

THE  EQUIVALENT  FREUNDLICH  PARAMETERS,  AND  THE  CONCENTRATION  RANGE 
WERE  THE  D-R  PREDICTION  AND  THE  FREUNDLICH  PREDICTION  ARE  WITHIN  10  % 


CHARACTER  NAME*25 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 
DOUBLE  PRECISION  N,NA,NB 
OPEN ( 7 , FILE* ' SPEQ . OUT ' ) 

READ  IN  THE  INPUT  DATA  FROM  THE  SCREEN 


WRITE (*,101) 

101 

FORMAT (2X, 'BEGIN  EXECUTION  OF 
write (*,5000) 

PROGRAM' 

.///) 

5000 

FORMAT (2X, 'ENTER  THE  COMPOUND 
READ(*,6500)  NAME 

NAME  ? 

') 

6500 

FORMAT (A2 5) 

WRITE (*,600) 
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600  FORMAT {2X, 'ENTER  ANTOINES  CONSTANT  A  ?  ') 

READ(*,*)  A 
WRITE (*,700) 

700  FORMAT (2X, 'ENTER  ANTOINES  CONSTANT  B  ?  ') 

READ(*,*)  B 
WRITE (*,800) 

800  FORMAT ( 2 X, 'ENTER  ANTOINES  CONSTANT  C  ?  ') 

READ(*,*)  C 
WRITE (*,900) 

900  FORMAT ( 2 X, 'ENTER  THE  CONSTANT  BETA  ?  ') 

READ(*,*)  BETA 
WRITE (*,100) 

100  FORMAT (2X, 'ENTER  THE  DENSITY  OF  THE  COMPOUND  IN  GM/CM**3  ?  ') 

READ(*,*)  DEN 
WRITE (*,110) 

110  FORMAT (2X, 'ENTER  THE  MOLECULAR  WEIGHT  OF  THE  COMPOUND  ?  ') 

READ(*,*)  WEIGHT 
WRITE (*,120) 

120  FORMAT (2X, 'ENTER  THE  SYSTEM  PRESSURE  IN  mm  HG  ?  ') 

READ(*,*)  PRESS 
WRITE (*,130) 

130  FORMAT (2X, 'ENTER  THE  SYSTEM  TEMPERATURE  IN  DEG.  C  ?  ') 

READ(*,*)  TEMPC 
WRITE (*,3434) 

3434  FORMAT ( 2 X, 'ENTER  THE  AVERAGE  INFLUENT  CONC.  IN  UG/L  ?  ') 

READ(*,*)  CONIN 
WRITE(*,1111) 

1111  FORMAT ( 2 X, 'ENTER  THE  INITIAL  CONC.  FOR  INTEGRATION  IN  UG/L  ?  ') 
READ(*,*)  CONC 

WRITE (*,1112) 

1112  FORMAT (2X, 'ENTER  THE  INTEGRATION  CONC.  STEP  IN  UG/L  ?  ') 
READ(*,*)  H 

C 

C  CALCULATE  THE  COMPOUNDS  VAPOR  PRESSURE  AT  THE  SYSTEMS  PRESSURE 
C  USING  THE  ANTOINES  EQUATION 
C 

TEMPK  »  TEMPC  +  273.0 

VP  »  10.0** (A  -  B/(TEMPC+C)) 

C  PRINT  OUT  THE  DATA  READ  INTO  THE  PROGRAM 
C 

WRITE(7,10)  NAME, A,  B,C,  BETA,  DEN,  WEIGHT,  PRESS,  TEMPC, CONIN 
10  FORMAT ('  COMPOUND  -  ',  A25,// 

$  '  THE  ANTOINES  EQTN.  A  IS  ',  F7.2,  / 

$  '  THE  ANTOINES  EQTN.  B  IS  ',  F7.2,  / 

$  '  THE  ANTOINES  EQTN.  C  IS  ',  F7.2,/ 

$  '  THE  VALUE  OF  BETA  IS  ',  F7.5,/ 

$  '  THE  LIQUID  DENSITY  IN  G/CC  IS  ',  F7.5,/ 

$  '  THE  MOLECULAR  WEIGHT  IS  ',  F6.2,/ 

$  '  THE  SYSTEM  PRESSURE  IN  MM  HG  IS  ',  F6.2,  / 

$  '  THE  SYSTEM  TEMPERATURE  IN  DEG  C  IS  ',  F6.2,/ 

$  '  INLET  GAS  CONCENTRATION  IN  ug/1  IS  ',  E12.5,//) 

WRITE (7, 20) 

20  FORMAT(2X,'  CONC  ug/1  ',3X,  'Qo  umol/gram' , 3X, '  N  ', 

$  3X, 'SPREADING  PRESSURE' ,/) 

CSC  INTIALIZE  THE  VALUES  FOR  THE  NUMERICAL  INTEGRATION 
C 

NA  -  O.ODO 
QA  -  O.ODO 
SUM  -  O.ODO 
EXC  «  l.OD-4 
I  «  1 
MAX  -  10 


55 


30  IF  (CONC  .LE.  CONIN)  THEN 
C 

C  CALCULATE  THE  SURFACE  LOADING  AT  THE  GAS  CONCENTRATION 
C 

PP  »  (CONC/WEIGHT  *  22.40D0  /  l.OD+06)  *  PRESS 
QB  s  0.46D-*-6*DEN/WEIGHT  * 

$  EXP(-1.33D-7*{TEMPK/BETA  *  DLOG(VP/PPJ ) **2 ) 

C 

C  CALCULATE  THE  VALUE  OP  din  C/d  In  q  AT  THE  GAS  CONCENTRATION 
C  THIS  IS  EQUAL  TO  THE  FREUNDLICH  N 
C 

NB  =  l.ODO  /  (2.66D-7  *  DLOG{VP/PP)  *  (TEMPK/BETA) **2 ) 

C 

C  TRAPEZOID  RULE 
C 

SUM  »  SUM  +  ((QB  -  QA)/2.0D0)  *  (NB  +  NA) 

NA  s  NB 
QA  -  QB 
C 

C  INCREMENT  THE  GAS  CONCENTRATION 
C 

IF  (I  .LE.  MAX)  THEN 
I  -  I  +  1 
CONC  »  CONC  +  H 
GO  TO  30 
ELSE 
C 

C  PRINT  OUT  INTERMEDIATE  VALUES  OF  SPREADING  PRESSURE 
C 

WRITE (7,40) CONC, QA, NA, SUM 
MAX  «  MAX  MAX 
IP  (CONC  .GT.  EXC)  THEN 
H  *  H  *  lO.ODO 
EXC  »  EXC  *  lO.ODO 

en'dif 

CONC  »  CONC  +  H 
PRINT*,  'CONC  IS  ',  CONC,QA 
I  »  I  +  1 
GO  TO  30 
END  IF 
ELSE 

WRITE (7,40) CONC, QA, NA, SUM 

40  FORMAT(2X,E12.5,3X,E12.5,3X,E12.5,3X,E12.5) 

C 

C  CALCULATE  THE  FREUNDLICH  PARAMETERS  AND  PRINT  OUT  THE  RESULTS 
C 

N  -  QA/SUM 

XK  “  QA/  ((CONC/WEIGHT)  **  N) 

WRITE(7,50)  SUM,N,XK 
50  FORMAT('  ' , // 

$  '  THE  SPREADING  PRESSURE  IS  *,E12.5,/ 

$  '  THE  VALUE  OF  1/N  IS  ',E12.5,/ 

$  '  THE  VALUE  OF  XK  IS  ',E12.5,//) 

ENDIF 

C 

C  CALCULATE  AND  PRINT  OUT  THE  LOWER  BOUNDS  FOR  THE  REGION  WHERE  THE 
C  ERROR  IS  LESS  THAN  10  %  BETWEEN  THE  D-R  PREDICTION  AND  THE 
C  FREUNDLICH  PREDICTION 
C 

CONC  -  CONIN 
QB  «  QA 
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o  o  o  o  o 


60 


IF(ABS((QB  -QA)/QB)  .LE.  0.10)  THEN 

PP  »  (CONC/WEIGHT  *  22.40D0  /  l.OD+6) ‘PRESS 
QB  »  0.46D+6*DBN/WEIGHT  * 

$  EXP(-1.33D-7*(TEMPK/BETA  *  DLOG(VP/PP) ) **2 ) 

QA  «  XK  *  (CONC/WEIGHT) **N 
CONC  »  CONC  -  0.001 
GO  TO  60 
ELSE 

WRITE(7,65)  CONC,QB,QA 
65  FORMAT ('  LOWER  CONCENTRATION  BOUNDS  IS  E12.5,'  ug/L',/ 

$  '  D-R  PREDICTED  LOADING  IN  umol/g  IS  E12.5,/ 

$  '  FREUNDLICH  LOADING  IN  lunol/g  IS  E12.5,/) 

ENDIF 

CALCULATE  AND  PRINT  OUT  THE  UPPER  BOUNDS  FOR  THE  REGION  WHERE  THE 
ERROR  BETWEEN  THE  D-R  PREDICTION  AND  FREUNDLICH  PREDICTION  IS 
LESS  THEN  10  % 


CONC  -  CONIN 
QB  -  QA 

70  IF(ABS((QB  -QA)/QB)  .LE.  0.10)  THEN 

PP  *  (CONC/WEIGHT  *  22.40D0  /  l.OD+6) ‘PRESS 
QB  =  0.46D+6‘DEN/WEIGHT  ‘ 

$  EXP(-1.33D-7‘(TEMPK/BETA  ‘  DLOG(VP/PP) )  “2 ) 

QA  »  XK  ‘  (CONC/WEIGHT)  “N 
CONC  *  CONC  +  0.001 
GO  TO  70 
ELSE 

WRITE(7,75)  CONC,QB,QA 

75  FORMAT ('  UPPER  CONCENTRATION  BOUNDS  IS  E12.5,  '  ug/L',/ 

$  '  D-R  PREDICTED  LOADING  IN  umol/g  IS  E12.5,/ 

$  '  FREUNDLICH  LOADING  IN  umol/g  IS  E12.5,/) 

ENDIF 

END 
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EXHIBIT  A-2 


COMPOUND  -  m-xylene 

THE  ANTOINES  EQTN.  A  IS  7.01 

THE  ANTOINES  EQTN.  B  IS  1462.00 

THE  ANTOINES  EQTN.  C  IS  215.00 

THE  VALUE  OF  BETA  IS  1.15830 

THE  LIQUID  DENSITY  IN  G/CC  IS  0.86400 

THE  MOLECULAR  WEIGHT  IS  106.00 

THE  SYSTEM  PRESSURE  IN  MM  HG  IS  760.00 

THE  SYSTEM  TEMPERATURE  IN  DEG  C  IS  25.00 

INLET  GAS  CONCENTRATION  IN  ug/1  IS  0.20000E-f03 


SPREADING  PRESSURE 


CONC  ug/1 

O.llOOOE+03 

0.21000E-f03 


Qo  umol/gram 

0.26874E-«-04 

0.26874E-«-04 


N 

0.92344E+01 

0.92344E+01 


0.12136E+05 

0.12136E-t-05 


THE  SPREADING  PRESSURE  IS  0 . 12136E-t-05 
THE  VALUE  OF  1/N  IS  0.22145E+00 
THE  VALUE  OF  XK  IS  0. 23098E-«-04 


LOWER  CONCENTRATION  BOUNDS  IS  0.15240E+03  ug/L 
D-R  PREDICTED  LOADING  IN  umol/g  IS  0.27814E+04 
FREUNDLICH  LOADING  IN  umol/g  IS  0.25032E-t-04 

UPPER  CONCENTRATION  BOUNDS  IS  0.69497E+03  ug/L 
D-R  PREDICTED  LOADING  IN  umol/g  IS  0.31845E+04 
FREUNDLICH  LOADING  IN  umol/g  IS  0.35029E^04 
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APPENDIX  B 


EVALUATION  OF  MULTICX)MPONENT  CXDMPETmVE  LOADING  IN  FIXED  BEDS 

The  program  ETMOD.FOR  is  a  modification  of  the  program  ET  furnished  by  J.  C.  Crittenden 
of  Michigan  Technology  University.  The  modification  enables  use  of  the  current  International 
Machine  and  Statistics  Library  (IMSL)  for  the  solution  of  a  set  of  mathematical  equations.  The 
current  IMSL  routine  used  in  this  calculation  is  DNEQNF.  The  use  of  this  program  is  described  in 
Crittenden  et  al.  (Reference  25).  The  necessary  input  data  are  described  below: 

Line  1 

Number  of  components  (N),  bed  void  fraction,  bulk  density  of  adsorbent,  flow  rate  of  fluid 
(gallons  per  minute  per  square  foot) 

Line  2  through  Line  N+l 

Component  name,  Freundlich  K  (micromole/gram)  (liters/micromole)’^,  Freundlich  1/n,  initial 
concentration  (micrograms/liter),  molecular  weight 

Typical  input  is  provided  in  Exhibit  B-1;  typical  output  is  found  in  Exhibit  B-2.  The  computer 
program  ETMOD.FOR  is  shown  in  Exhibit  B-3.  Note  that  the  input  flow  rate  is  in  gallons  per 
minute  per  square  foot.  The  program  was  originally  written  for  water  treatment.  However,  the 
program  may  be  used  for  air  treatment  by  (1)  converting  typical  air  flow  rate  units  (e.g.,  cubic  feet 
per  minute)  to  those  specified  for  the  input  flow  rate,  and  (2)  inputting  Freundlich  parameters  that 
apply  for  gas-phase  adsorption. 
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EXHffiITB-1 


TYPICAL  INPUT 

NOTE:  This  is  liquid-phase  data;  however,  gas-phase  data  can  be  used. 
SAMPLE  INPUT  -  ETJDAT 


7,0.431,0.457,1.931 

’M-XYLENEM044.0, 0.2458,5.2, 106.17 

’OP-XYLENE’,  895.0,0.2587,9.3,106.17 

’ETHYL  BENZENE’,  714.4,0.2953,4.5,106.17 

TOLUENE’, 475.0,0.3282, 19.3,92.15 

’TETRACHLOROETHENE’,435.0,0.3847,37.6,165.83 

’TRICHLOROETHENE’,192.0,0.4382,47.9,131.29 

’CIS-DICHLOROETHENE’,46.9,0.5562,70.9,96.93 
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EXHIBIT  B-2 


TYPICAI.  OUTPUT 


EQUILIBRIUM  THEORY  RESULTS 


NUMBER  OF  COMPONENTS  : 

5 

BED 

VOID 

FRACTION  : 

0.430 

BULK  DENSITY  OF  ADSORBENT 

(g/cm**3)  : 

0.480 

FLOWRATE 

(gpm/ft’**2)  : 

224.400 

BED  VOLUMES  FED 

VELOCITY  OF  WAVE 

TREATMENT  CAPACITY 

TO  BREAKTHROUGH 

(cm/sec) 

(mg  CARBON/L  WATER) 

ZONE( 

1) 

334614.5 

0.45541907E-04 

1.4345 

ZONE( 

2) 

312381.6 

0.48783232E-04 

1.5366 

ZONE( 

3) 

281513.9 

0.54132259E-04 

1.7051 

ZONE( 

4) 

250931.5 

0.60729666E>04 

1.9129 

ZONE( 

5) 

177584.1 

0.85812782E-04 

2.7029 

COMPONENT*  M, O , P-XYLENES 


FREUNDLICH  K  (um/g) * (L/um) **l/n  : 
FREUNDLICH  1/n  : 

INITIAL  CONCENTRATION  (ug/L)  : 
MOLECULAR  WEIGHT  : 

SINGLE  SOLUTE  TREATMENT  CAPACITY 
(mg  CARBON/L  WATER)  : 
DIMENSIONLESS  BED  LENGTH 
FOR  ZONE(  1)  t 


2390.00 

0.1920 

356.0000 

106.1700 

1.1122 

0.53071239 


C(ug/L) 

Q(ug/g) 

C/Co 

DG 

QAVE 

Cio  (ug/L) 

ZONE( 

1) 

356.000 

248172.46 

1.0000 

.778E+06 

248172.46 

0.521E+03 

ZONE( 

2) 

0.000 

0.00 

0.0000 

.OOOE+OO 

231683.03 

0.467E+03 

ZONE( 

3) 

0.000 

0.00 

0.0000 

. OOOE+00 

208789.49 

0.806E+02 

ZONE( 

4) 

0.000 

0.00 

0.0000 

.OOOE+00 

186107.51 

0.517E+01 

ZONE( 

5) 

0.000 

0.00 

0.0000 

.OOOE+00 

131708.20 

0.134E-02 

COMPONENT*  ETHYL  BENZENE 


FREUNDLICH  K  (um/g) * (L/uro) **l/n  ;  2481.00 

FREUNDLICH  1/n  :  0.2194 

INITIAL  CONCENTRATION  (ug/L)  ;  178.0000 

MOLECULAR  WEIGHT  :  106.1700 

SINGLE  SOLUTE  TREATMENT  CAPACITY 

(mg  CARBON/L  WATER)  :  0.6033 

DIMENSIONLESS  BED  LENGTH 

FOR  ZONE(  2)  :  0.03777206 


C(ug/L) 

Q(ug/g) 

C/Co 

DG 

QAVE 

Cio  (ug/L) 

ZONE( 

1) 

178.000 

97709.62 

1.0000 

.613B+06 

97709.62 

0.662E+03 

ZONE( 

2) 

569.459 

370601.29 

3.1992 

.726E+06 

115841.48 

0.602E+03 

ZONE( 

3) 

0.000 

0.00 

0.0000 

.OOOE'fOO 

104394.71 

0.129E+03 

ZONE( 

4) 

0.000 

0.00 

0.0000 

.OOOE+00 

93053.73 

0.117E+02 

ZONE( 

5) 

0.000 

0.00 

0.0000 

.OOOE+00 

65854.08 

0.847E-02 
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COMPONENT*  TOLUENE 


ZONE( 

ZONE( 

ZONE( 

ZONE( 

ZONE( 


FREUNDLICH  K  (uin/g)  *  (L/um)  **l/n  ; 
FREUNDLICH  1/n  : 

INITIAL  CONCENTRATION  (ug/L)  : 

molecular  weight  : 

SINGLE  SOLUTE  TREATMENT  CAPACITY 
(mg  CARBON/L  WATER)  : 
DIMENSIONLESS  BED  LENGTH 
FOR  ZONE(  3)  : 


1) 

C(ug/L) 

62.530 

Q(ug/g) 

10453.11 

C/Co 

1.0000 

2) 

65.897 

12800.62 

1.0539 

3) 

498.542 

294363.26 

7.5654 

4) 

0.000 

0.00 

0.0000 

5) 

0.000 

0.00 

0.0000 

2029.00 

0.2779 

62.5300 

92.1400 

0.3725 

0.06233369 

DG  QAVE  Cio  (ug/L) 

.187E+06  10453.11  0.189E+04 

.217E+06  10609.09  0.175E+04 

.659E+06  38647.99  0.520E+03 

.OOOE+00  34449.44  0.779E+02 

.OOOE+00  24379.85  0.259E+00 


COMPONENT*  BENZENE 

FREUNDLICH  K  (um/g) * (L/um) **l/n  : 

FREUNDLICH  1/n  : 

INITIAL  CONCENTRATION  (ug/L)  : 

MOLECULAR  WEIGHT  : 

SINGLE  SOLUTE  TREATMENT  CAPACITY 
(mg  CARBON/L  WATER)  : 

DIMENSIONLESS  BED  LENGTH 
FOR  ZONE(  4)  : 

C/Co  DG  QAVE  Cio  (ug/L) 

1.0000  .864E+05  3555.72  0.345E+04 

1.0167  .978B+05  3591.28  0.329E+04 

1.1805  .194E+06  4182.55  0.150E+04 

7.8474  .596E+06  28817.53  0.438E+03 

0.0000  .OOOE+00  20394.16  0.108E+02 


ZONE(  1) 
ZONB(  2) 
ZONE (  3 ) 

ZONE (  4 ) 


C(ug/L) 

45.930 

46.698 

55.124 

432.583 


Q(ug/g) 

3555.72 

4090.92 

9574.93 
230949.58 

0.00 


1427.00 

0.4288 

45.9300 

78.1100 

0.5174 

0.07688140 


COMPONENT*  N-PENTANE 


FREUNDLICH  K  (um/g) * (L/um) **l/n  :  518.00 

FREUNDLICH  1/n  :  0.4876 

INITIAL  CONCENTRATION  (ug/L)  :  37.4400 

MOLECULAR  WEIGHT  ;  72.1500 

SINGLE  SOLUTE  TREATMENT  CAPACITY 

(mg  CARBON/L  WATER)  :  1.3794 

DIMENSIONLESS  BED  LENGTH 

FOR  ZONE(  5)  :  0.29230047 


ZONE( 

1) 

C(ug/L) 

37.440 

Q(ug/g) 

440.13 

C/Co 

1.0000 

DG 

.131E+05 

QAVE 

440.13 

Cio  (ug/L) 
0.210E+05 

ZONE( 

ZONE( 

ZONE( 

ZONE( 

2) 

37.520 

496.19 

1.0021 

.148E+05 

443.85 

0.201E+05 

3) 

38.156 

909.92 

1.0169 

.266E+05 

489.91 

O.lOlE+05 

4) 

5) 

41.025 

2592.67 

1.0752 

.705E+05 

718.34 

0.342E+04 

132.067 

50186.75 

3.2192 

.424E+06 

15177.98 

0.132E+03 
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EXHIBIT  B3 


C 

C 

C  ********EQUILIBRIUM  THEORY  PROGRAM******** 

C 

C 

C  THIS  PROGRAM  CALCULATES  MULTICOMPONENT  BREAKTHROUGH  FOR  FIXED 
C  BED  ADSORBERS.  THE  PROGRAM  ASSUMES  NO  MASS  TRANSFER  RESISTANCE. 

C  IDEAL  ADSORBED  SOLUTION  THEORY  IS  USED  TO  PREDICT  COMPETITION. 

C 

C 

C  PROGRAM  WRITTEN  BY:  THOMAS  FRANCIS  SPETH 

C  PROGRAM  ALTERED  BY  R  M  COUNCE  NOV  29,1989 

C  THANKS  TO:  PAUL  LUFT,  DAVID  HAND,  AMD  DR.  JOHN  CRITTENDEN 
C  FOR  THEIR  PREVIOUS  WORK. 

C 

C  VARIABLE  DEFINITIONS 

C 

C  BVF  >  BED  VOID  FRACTION 

C  C  x  LIQUID  PHASE  CONCENTRATION  (ug/L) 

C  CH  «  WORKING  CHARACTER 

C  CHAR  X  NAME  OF  THE  COMPONENTS  (TEN  LETTERS) 

C  CO  X  initial  concentrations  (ug/L) 

C  DEN  X  bulk  density  OF  ADSORBENT  (g/cm**3) 

C  DGX  X  DIMENSIONLESS  GROUP  X:  USED  TO  FIND  STRONGEST  COMPONENT 

C  DGY  -  DIMENSIONLESS  GROUP  Y:  USED  TO  FIND  STRONGEST  COMPONENT 

C  FCN  X  SUBROUTINE  THAT  SETS  UP  THE  NON-LINEAR  EQUATIONS 
C  FCS  X  C/CO 

C  FLRT  X  FLOW  RATE  (GPM/FT**2) 

C  FNORM  X  OUTPUT:  SUM  OF  THE  RESIDUALS 
C  I  X  COUNTER 

C  lAST  X  SUBROUTINE  TO  ACCOUNT  FOR  COMPETITIVE  EFFECTS 
C  lER  X  OUTPUT:  ERROR  PARAMETER 

C  ITMAX  X  MAXIMUM  NUMBER  OF  ITERATIONS 

C  IX  X  USED  TO  KEEP  TRACK  OF  STRONGEST  COMPONENT 
C  J  X  COUNTER 

C  K  X  COUNTER 

C  L  X  COUNTER  FOR  ERROR  FIXING 

C  M  X  COUNTER 

c  MW  X  molecular  weight 

c  N  X  number  of  components  total 

C  NN  X  NUMBER  OF  COMPONENTS  IN  A  ZONE 

c  NS  X  nsig  input 

C  NSIG  X  NUMBER  OF  DIGITS  OF  ACCURACY  DESIRED  IN  THE  COMPUTED  ROOT 

C  OATS  X  bed  volumes  FED 

C  PAR  X  parameter  SET 

C  PAR(1  to  N)x  FREUNDLICH  K  VALUES 
C  PAR(10  to  10-»-N)x  FRUENDLICH  N  VALUES 

C  PAR (20  to  20+N)-  INITIAL  CONCENTRATIONS 

C  PAR(30)x  VELOCITY  OF  THE  WAVE:  VW  (cm/s) 

C  PAR(35)x  VELOCITY  OF  FLOW:  VF  (cm/s) 

C  PAR (40  to  40+1 )x  CALCULATED  LIQUID  CONCENTRATIONS 

C  PAR (60  to  60+1 )«  Q's  OF  THE  PREVIOUS  WAVE 

C  PAR (80  to  80+1 )>  C'B  OF  THE  PREVIOUS  WAVE 


c 

Q 

SB 

SOLID  PHASE  CONCENTRATION  (ug/g) 

c 

QAVE 

m 

AVERAGE  Q  IN  ZONE 

c 

SSTC 

m 

SINGLE  SOLUTE  TREATMENT  CAPACITY  (mg  C/L  WATER) 

c 

SUM 

m 

USED  TO  CALCULATE  VW  AND  OATS 

c 

VF 

m 

VELOCITY  OF  FLOW  (cm/s) 

c 

VW 

m 

VELOCITY  OF  WAVE  (cm/s) 

c 

WK 

m 

WORK  VECTOR:  LENGTHxN*(3*N+15)/2 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


10 

c 

c 

c 


c 

c 

c 


20 

c 

c 

c 


c 


X  =  ONE  DIMENSIONAL  SOLID-PHASE  CONCENTRATION  (um/g) 

XK  »  FRUENDLICH  K's  (um/g) * ( (L/um) **1/N 

XN  »  FRUENDLICH  1/n  'a 

ZSQ  >  COMMON  BLOCK 

ZZ  «  VARIABLE  USED  TO  CALCULATE  INITIAL  Q'a 

ZZZ  >  DIMENSIONLESS  BED  LENGTH 


SAMPLE  INPUT 
N,BVF,DEN,FLRT 

CHAR(I) ,XK(I) ,XN(I),CO(I),MW(I) 
•  •••  ••  ••  ••  •• 
CHAR( J) ,XK( J) ,XN( J) ,CO( J) ,MW( J) 
NS 


DIMENSIONS 

IMPLICIT  DOUBLE  PRECISION  (A-H,  O-Z) 

CHARACTER*40  CHAR(20),CH 
DOUBLE  PRECISION  MW (20) 

DIMENSION  XN(20) ,XK(20) ,CO(20,20),C(20,20) ,Q(20,20) ,ZZZ(20) ,VW(20) 
$,DGy(20,20) , OATS (20) ,QAVE(20,20) ,SSTC(20) ,WK(200) ,X(20) 
$,XGUESS(20) ,FCS(20,20) ,CIO(20,20) 

EXTERNAL  FCN 

COMMON  /ZSQ/  BVF,DEN,M,PAR(100) 

OPEN  FILES 

OPEN ( 4 , FILE- ' ET . DAT ' , STATUS- ' OLD ' ) 

OPEN ( 7 , FILE- ' ET . OUT ' , STATUS- ' NEW ' ) 

READ  IN  DATA 

PRINT* ,  ' READ ING  DATA ' 

REA0(4,*)  N,BVF,DBNrFLRT 
DO  10  I»1,N 

READ (4,*)  CHAR(I) ,XK(I) ,XN(I) ,CO(I,l) ,MW(I) 

CO(I,l)-CO(I,l)/MW(I) 

XN(I)-1.0D0/XN(I) 

CONTINUE 


CHANGE  UNITS 

VF«FLRT*0.067910D0/BVF 
DEN-DEN* 1000 . 000 

SET  ZONE  ONE  CONCENTRATIONS  TO  ZERO 

DO  20  I-1,N 
VW(I)-0.0D0 
PAR(60-«>I)-0.0D0 
PAR ( 80+1 )»0.0D0 
CONTINUE 


SOLVE  FOR  EACH  ZONE  SEPARATELY 


DO  100  J»1,N 

PRINT*,  'SOLVING  FOR  ZONE  ',J 
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L»0 

M=J 

NS>9 

NN»N+1-J 

ZZ-l.ODO 

NSIG^NS 

SUM^O.ODO 

IRRR£L>1.0E-05 

C 

C  CALCULATE  INITIAL  GUESSES  OF  Q'S 

C 

21  DO  22  I»1,N 

Q(I,J)»ZZ*XK(I)*CO(I, J)**(1.0D0/XN(I) ) 

22  CONTINUE 
C 

C  PUT  Q  INTO  ONE-DIMENSIONAL  FORM 

C 

DO  24  I-1,N 
XGUESS(I)«Q(I, J) 

24  CONTINUE 

C 

C  SET  lAST  PARAMETERS 

C 

DO  26  I>1,NN 
XGUESS(I)«X(M-1+I) 

PAR(I)»XK(M-1+I) 

PAR( lO+I ) »XN(M-1+I ) 

PAR(20+I)»CO(M-1+I, J) 

PAR(60+I)-PAR(60+M-1+I) 

PAR  ( 80+ 1 )  =*PAR  ( 80+M- 1+ 1 ) 

26  CONTINUE 

PAR<30)»VW(J-1) 

PAR(35)-VF 

ITMAX-100 

C 

CALL  DNEQNF ( FCN , IRRREL , NN , ITMAX, XGUESS , X , FNORM ) 

C 

C  FIX  ANY  ERRORS 

C 

IF  (lERCDO  .EQ.  1  .OR.  IERCO()  .EQ.  3)  THEN 
IF  (L  .EQ.  0)  THEN 
ZZ>2.0D0*ZZ 
L«L+1 
GOTO  21 
ENDIF 

IF  (L  .EQ.  1)  THEN 
ZZ-3.0D0+ZZ 
L«L+1 
GOTO  21 
ENDIF 

IF  (L  .EQ.  2)  THEN 
ZZ-ZZ/20.0D0 
L-L+1 
GOTO  21 
ENDIF 

IF  (L  .EQ.  3)  THEN 
WRITE(7,28)  J 

28  FORMAT( IX, 'THERE  IS  A  PROBLEM  WITH  THE  INITIAL  CONCENTRATIONS', 

$/,'THAT  THE  PROGRAMS  FIXING  ROUTINE  DID  NOT  HELP.  ZONE=',I2) 
ENDIF 
ENDIF 
C 
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IF  (lERCOO  .EQ.  2)  THEN 
NSIG-NSIG-1 
IF  (NSIG  .LT.  O)  THEN 
WRITE (7, 29)  J 

29  FORMAT(lX, 'THE  NimSER  OF  SIGNIFICANT  FIGURES  HAS  DROPPED  BELOW 
$ZERO.  THERE  ARE  NO  RESULTS  FOR  ZONE'/ 12) 

ENDIF 
GOTO  21 
ENDIF 
C 

C  SET  X  TO  TWO-DIMENSIONAL  OUTPUT  FOR  PRINT  OUT 

C 

IF  (M  .GT.  1)  THEN 
DO  30  I»1,M-1 
Q(I, J)»0.0D0 

30  CONTINUE 
ENDIF 

DO  31  I«1,NN 
Q(1+M-1,J)«X(I) 

31  CONTINUE 
C 

C  CALCULATE  THE  LIQUID  CONCENTRATIONS 

C 

IF  (M  .GT.  1)  THEN 
DO  33  1-1, M-1 
C(I,J)-0.0D0 

33  CONTINUE 
ENDIF 

DO  34  I-1,NN 
C(1+M-1,J) -PAR (40+1) 

34  CONTINUE 
C 

C  DETERMINE  THE  STRONGEST  COMPONENT  IN  ZONE  J 

C 

DGX-O.ODO 
DO  35  I-M,N 

DG-DEN*Q(I, J)/(C(I,J)*BVF) 

IF  (DG  .GT.  DGX)  THEN 
DGX-DG 
IX-I 
ENDIF 

35  CONTINUE 
C 

C  SET  STRONGEST  COMPONENT  TO  ZONE  J 

C 

CH-CHAR(IX) 

CHAR(IX)-CHAR(J) 

CHAR(J)-CH 

C 

WM-MW(IX) 

MW(IX)=MW(J) 

MW( J)-WM 
C 

XXK-XK(IX) 

XK(IX)-XK(J) 

XK(J)-XXK 

C 

XXN-XN(IX) 

XN{IX)-XN(J) 

XN(J)-XXN 


C 
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DO  37  K=1,J 
XCO=CO(IX,K) 

CO(IX,K)=CO( J,K) 

CO(J,K)*XCO 

C 

XC*C(IX,K) 

C(IX,K)=C{ J,K) 

C{J,K)=XC 

C 

XQ=Q(IX,K) 

Q(IX,K)*Q(J,K) 

Q{J,K)=XQ 

37  CONTINUE 
C 

C  SET  C'a  AND  Q's  FOR  NEXT  ZONE 

C 

DO  38  1=1, N 
PAR(60+I)=Q(I, J) 

PAR(80+I)=C(I, J) 

38  CONTINUE 
C 

C  CALCULATE  VELOCITY  OF  THE  WAVE  FOR  ZONE  J 

C 

IF  (J  .EQ.  1)  THEN 

VW(J)=VF*BVF*CO(l, J)/(Q(1, J)*DEN+C(1, J)*BVF) 

ENDIF 

IF  (J  .GE.  2)  THEN 
SUM*(Q(J,1)*DEN+BVF*C( J,l) )*VW(1) 

ENDIF 

IF  (J  .GT.  2)  THEN 
DO  40  K=2,J-1 

SUM»SUM+( (Q(J,K)*DEN+BVF*C(J,K) )*(VW(K)-VW(K-1) ) ) 

40  CONTINUE 

ENDIF 

IF  (J  .GE.  2)  THEN 

VW( J)-(BVF*VF*CO( J, J)-SUM+(Q(J, J)*DEN+BVF*C( J, J) )*VW( J-1) 
$)/(Q(J» J)*DEN+BVF*C(J, J) ) 

ENDIF 

C 

C  SET  Co  FOR  NEXT  ZONE 

C 

DO  50  I=J+1,N 
CO(I,J+l)»C(I, J) 

50  CONTINUE 

DO  60  1=1, J 
CO(I, J+1)=0.0D0 
60  CONTINUE 

C 

100  CONTINUE 
C 

C  CALCULATE  BED  VOLUMES  FED 

C 

PRINT*,  'CALCULATING  VARIOUS  PARAMETERS' 

DO  110  1*1, N 

SUM=(Q(I,1)*DEN+C(I, J)*BVF)*VW(1) 

IF  (I  .GE.  2)  THEN 
DO  105  K=2,I 

SUM*SUM+(Q(I,K)*DEN+C( J,K)*BVF)*{VW(K)-VW(K-1) ) 

105  CONTINUE 

ENDIF 

OATS(I)-SUM/(CO(I,I)*VW(I) ) 

110  CONTINUE 
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c 

C  CALCULATE  Q  TOTAL  AVERAGE  FOR  EACH  ZONE 

C 

DO  120  Ial,N 
QAVE(I,1)=Q(I,1) 

120  CONTINUE 

DO  140  1=1, N 
DO  130  J=2,N 

QAVE(I, J)=(QAVE(I, J-1)*VW(J-1)+Q(I, J)*(VW( J)-VW(J-l) ) )/VW(J) 
130  CONTINUE 
140  CONTINUE 
C 

C  CALCULATE  DIMENSIONLESS  BED  LENGTH 

C 

DO  150  1=1, N 
IF  (I  .EQ.  1)  THEN 
ZZZ(1)=0.0D0 
ELSE 

ZZZ(I)=VW(I-1)/VW(N) 

ENDIF 

ZZZ(I)=VW(I)/VW(N)-ZZZ(I) 

150  CONTINUE 
C 

C  CALCULATE  SINGLE  SOLUTE  TREATMENT  CAPACITY 

C 

DO  160  I-1,N 

SSTC ( I ) =1000 . 0D0*CO ( I , 1 ) ** ( 1 . ODO-1 . ODO/XN ( I ) ) /XK ( I ) 

160  CONTINUE 
C 

C  CALCULATE  DG'8  FOR  PRINT  OUT 

C 

DO  180  J=1,N 
DO  170  1=1, N 

DGY(I, J)=DEN*Q(I, J)/(C(I,J)*BVP) 

IP  (C(I,J)  .Eg.  0.0)  DGY(I,J)=0.0 
170  CONTINUE 

180  CONTINUE 
C 

C  CALCULATE  C/Co  FOR  PRINT  OUT 

C 

DO  195  J=1,N 
DO  190  1*1, N 
FCS(I, J)=C(I, J)/CO(I, J) 

IF  (CO(I,J)  .EQ.  0.0)  FCS(I,J)=0.0 
190  CONTINUE 

195  CONTINUE 
C 

C  CALCULATE  Cio 

C 

DO  250  J»1,N 
CISUM=0.0D0 
DO  230  1=1, N 
C1SUM»CISUM+XN ( I ) *Q ( I , J ) 

230  CONTINUE 

DO  240  1=1, N 

C10(I, J)=(CISUM/XN(I)/XK(I) )**XN(I) 

240  CONTINUE 

250  CONTINUE 
C 
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PUT  INTO  ug/L  UNITS 


C 
C 

DO  210  J=1,N 
DO  200  1^1, N 
CO(I, J)=CO(I, J)*MW(I) 

CIO(I, J)=CIO(I, J)*MW(I) 

Q(I,J)=Q(I,J)*MW(I) 

QAVE  ( I ,  J )  =Q  AVE  ( I ,  J )  *  MW  ( I ) 

C(I,J)=C(I, J)*MW(I) 

200  CONTINUE 
210  CONTINUE 
C 

C  PRINT  RESULTS 

C 

PRINT*,  'PRINTING  OUTPUT' 

WRITE (7, 900) 

900  FORMAT(//, IX, TIO, 'EQUILIBRIUM  THEORY  RESULTS') 

WRITE (7, 910)  M,BVF, DEN/1000. 0D0,FLRT 
910  FORMAT( IX, /A TIO, 'NUMBER  OF  COMPONENTS  :',T54,I2,/, 

$T10,'BED  VOID  FRACTION  : ' , T55 , F5 . 3 , / , 

$T10,'BULK  DENSITY  OF  ADSORBENT  (g/cm**3)  : ' ,T50, FIO. 3 , / , 

$T10, 'FLOWRATE  (gpfn/ft**2)  : '  ,T50, FIO. 3, // ) 

WRITE (7, 920) 

920  FORMAT (T19, 'BED  VOLUMES  FED' ,T36, 'VELOCITY  OF  WAVE' ,T54, 

$' TREATMENT  CAPACITY',/, 

$T19,'TO  BREAKTHROUGH' ,T40, ' (cm/sec) ',T54, ' (mg  CARBON/L  WATER)',/) 
DO  940  I=-1,N 

WRITE(7,930)  I,OATS(I) ,VW(I) ,DEN*1000.0D0/OATS{I) 

930  FORMAT(T8, ' ZONE ( ' , T13 , 12 , T15 , ' ) ' ,T17, F13 . 1 , T36, E15 . 8, T52 , F15 . 4 ) 

940  CONTINUE 

DO  1000  1=1, N 

WRITE (7, 950)  CHAR(I) ,XK(I) ,1.0D0/XN(I) ,CO(l, 1) ,MW(1) ,SSTC(I) , I, 

$  ZZZ(I) 

950  FORMAT(//, 'COMPONENT-  ',A40,//, 

S  TIO, 'FREUNDLICH  K  (um/g) * (L/um) **l/n  : ' ,T45,F10.2, /, 

$  TIO, 'FREUNDLICH  1/n  : ' ,T45 ,F10. 4, / , 

$  TIO, 'INITIAL  CONCENTRATION  (ug/L)  : ' , T45 , FIO . 4 , / , 

$  TIO, 'MOLECULAR  WEIGHT  : ' ,T45,F10.4, /, 

$  TIO, 'SINGLE  SOLUTE  TREATMENT  CAPACITY',/, 

$  T12,'(mg  CARBON/L  WATER)  ; ' ,T40, F15 . 4, / , 

$  TIO, 'DIMENSIONLESS  BED  LENGTH', ^ 

$  T12,'FOR  ZONE(',T21,I2,T23,')  ; ' ,T45,F10.8, //, 

$  T12, 'C(ug/L) ' ,T24, 'Q(ug/g) ' ,T35, 'C/Co' ,T46, 'DG' ,T56, 'QAVE' ,T63, 

$  'Cio  (ug/L)') 

DO  970  J-1,N 

WRITE (7, 960)  J,C(I, J) ,Q( I , J) , FCS ( I, J) ,DGY(I, J) ,QAVE(I, J) , 

$  CIO(I,J) 

960  FORMAT ( ' ZONE ( ' , T7 , 12 , T9 , ' ) ' , Til , F9 . 3 , T22 , FIO . 2 , T33 , 

SF7 . 4 , T43 , E8 . 3 , T52 , FIO . 2 , T63 , ElO . 3 ) 

970  CONTINUE 
1000  CONTINUE 
STOP 
END 
C 
C 

C  SUBROUTINE  FCN 

C 

C  THIS  SUBROUTINE  WILL  SET  UP  THE  EQUATIONS  THAT  WILL  BE  USED  IN 
C  THE  ZSPOW  SUBROUTINE. 

C 

C 
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SUBROUTINE  FCN{X,F,NN) 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-2) 

DIMENSION  X(NN),F(NN) 

COMMON  /ZSQ/  BVF,DEN,M,PAR(100) 

PRINT*,  'M*',M 
QT=0 . ODO 
QNQ=0 . ODO 
DO  1010  1*1, NN 
QT=QT+X(I) 

QNQ*QNQ+PAR ( 10+1 ) *X ( I ) 

1010  CONTINUE 
C 

C  CALCULATE  F(I) 

C 

IF  (M  .EQ.  1)  THEN 

DO  1020  I»1,NN  ^ 

F ( I ) — PAR( 20+1 ) +X ( I ) /QT* (QNQ/PAR( 10+1 ) /PAR ( I ) ) *  *PAR ( lO+I ) 

1020  CONTINUE 
ENDIF 

IF  (M  .GT.  1)  THEN 
DO  1030  1=1, NN 

F(I)— X{I)/QT*(QNQ/PAR(10+I)/PAR(I) )**PAR(10+I)+( (X(I)-PAR(60+I) 
$)*DEN*PAR{30) )/( (PAR(35)-PAR(30) )*BVF)+PAR(80+I) 

1030  CONTINUE 
ENDIF 
C 

C  CALCULATE  LIQUID  CONCENTRATION 

C 

IF  (M  .EQ.  1)  THEN 
DO  1040  1=1, NN 

PAR(40+I)=X(I)*(QNQ/(PAR(10+I)*PAR(I)) )**PAR(10+I)/QT 

PRINT*,  'X(I)  ',X(I),'  F(I)  ',F(I) 

1040  CONTINUE 
ENDIF 

IF  (M  .GT.  1)  THEN 
DO  1050  1=1, NN 

PAR(40+I)=( (X(I)-PAR(60+I) ) *DEN*PAR( 30 ) ) / ( (PAR(35)-PAR(30) ) 
$*BVF)+PAR(80+I) 

PRINT*,  'X(I)  ',X(I),'  F(I)  ',F(I) 

1050  CONTINUE 
ENDIF 
RETURN 
END 
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APPENDIX  C 


SPREADSHEET  VARIABLES 

The  following  notation  table  applies  specifically  to  the  terminology  used  in  the  spreadsheets  for 
this  project.  It  is  important  to  recognize  that  the  spreadsheet  notation  often  is  not  the  same  as  the 
notation  used  in  the  preceding  sections  of  text.  This  arises  because  of  typical  software  restrictions 
as  to  acceptable  notation.  Because  of  this,  all  of  the  following  spreadsheet  notation  is  clearly  and 
explicitly  defined. 
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TABLE  C-1.  NOTATION  TABLE  -  DESIGN  AND  CXDST  ESTIMATING 


NAMES 

[bl 

ITEM  DESCRIPTION  [c] 

GASP 

N 

Catalytic  Unit  Size  Safety  Factor 

CD* 

N 

Carbon  Bulk  Density,  pounds/ft^ 

CP 

N 

Packing  Factor  (Treybal,  3rd  Ed.) 

CRT 

N 

Carbon  Recycle  Interval,  days 

CSP 

N 

Carbon  Use  Safety  Factor 

CUP 

N 

Contaminant  Cleanup  Factor 

EC* 

N 

Electricity  Cost,  $/kWh 

EFFA* 

N 

Electrical  Efficiency,  Air  Blowers 

EFFH* 

N 

Electrical  Efficiency,  Air  Heaters 

EFFW* 

N 

Electrical  Efficiency,  Water  Pumps 

EIFC* 

N 

Equip.  Inst  Cost  Factor,  Carbon  &  Catalysis 

EIFS* 

N 

Equipment  Installation  Cost  Factor,  Stripper 

P 

N 

Fanning  Friction  Factor  for  Water  Pipe 

FAP 

N 

Freundlich  Adjustment  Factor  (BENZENE  Only) 

FC* 

N 

Fuel  Oil  Costs,  S/Vf  BTU 

FP 

N 

Flooding  Factor  (SEE  ACCOMPANYING  TEXT) 

HGW* 

N 

Groundwater  Depth,  feet 

HPB* 

N 

Height  of  Each  Packed  Bed,  Feet 

HPL* 

N 

Horizontal  Pipe  Length,  feet 

I* 

N 

Average  Annual  Inflation  Rate 

J* 

N 

Correlation  Factor  Term  (Treybal,  3rd  Ed.) 

KP 

N 

Consolidated  Friction  Loss  Coefficient 

LP 

N 

Operating  Cycle  Load  Factor 

MMP 

N 

Material  (SS)  Factor,  Stripper 

OHR* 

N 

Overhead  Rate,  on  Expense 
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TABLE  C-1.  NOTATION  TABLE  -  DESIGN  AND  COST  ESTIMATING  (CONTINUED) 


NAMES 

Ibl 

ITEM  DESCRIPTION  [c] 

OM* 

N 

Operations  &  Maintenance  Factor 

PAT* 

N 

Gas  Stream  Pressure,  Atmospheres 

PPBP 

N 

Inlet  Contaminant  Concentration,  ppb 

PPBO 

N 

Outlet  Contaminant  Concentration,  ppb 

S* 

N 

Stripping  Factor  (Greater  Than  Unity) 

SP 

N 

Height  Design  Safety  Factor 

SIZE* 

B 

Size  of  Stripper  Packing,  mm  (inches) 

TP 

N 

Fahrenheit  Temperature 

TYPE* 

A 

Name/Type  of  Stripper  Packing 

VAR* 

N 

Independent  Variable,  Current  Value 

VI* 

A 

Name  of  Independent  Variable 

VIN* 

N 

Initial  Value  of  Independent  Variable 

VNUM* 

N 

Number  of  Parametric  Steps  to  Run 

VSTEP 

N 

Size  of  Each  Independent  Variable  Step 

WGM* 

N 

Water  Feed  Rate,  gallons/minute 

ABS 

N 

Abscissa  for  Pressure  Gradient  (Treybal,  3rd  Ed.) 

ABS2 

N 

Treybal  Abscissa,  Table  5 

ABSB 

N 

ABS  Individual  for  Benzene 

ABST 

N 

ABS  Individual  for  Trichloroethylene 

ACAT 

N 

Air  Flow  Rate,  ftVminute,  BENZENE 

ACFM 

N 

Air  Flow  Rate,  ftVminute 

ACFMl 

N 

ACFM  Individual  for  Benzene 

ACFM2 

N 

ACFM  Individual  for  Trichloroethylene 

ACFMB 

N 

Copied  from  ACFMl 

AD 

N 

Air  Density,  pounds/ft^ 
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TABLE  C-1.  NOTATION  TABLE  -  DESIGN  AND  CXDST  ESTTMATING  (CONTINUED) 


NAMES 

[b] 

ITEM  DESCRIPTION  [c] 

ALR 

N 

Air  Loading  Rate  at  Flooding,  (pounds/minute)/ft“ 

AMCF 

N 

Air,  pound  raoles/ft^ 

AORD 

N 

Adjusted  Ordinate,  FF  *  ORD 

APM 

N 

Air  Flow  Rate,  pounds/minute 

APMl 

N 

APM  Individual  for  Benzene 

APM2 

N 

APM  Individual  for  Trichloroethylene 

APMB 

N 

Air  Flow  Rate  BENZENE,  pounds/minute 

APMM 

N 

Air  Flow  Rate,  pound  moles/minute 

AREA 

N 

Stripper  Cross  Section  Area,  ft^ 

AVI 

N 

Air  Viscosity,  Centipoise 

AV2 

N 

Air  Viscosity,  pounds/feet  second 

BC 

N 

Carbon  Breakthrough  Capacity,  pounds/pounds 

CAPCATB 

N 

Catalysis  Unit  Installed  Capital  Cbst,  BENZENE 

CAPCB 

N 

Carbon  Unit  Installed  Capital  Costs,  BENZENE 

CAPRB 

N 

Centrifugal  Contactor  Installed  Capital  Costs,  BENZENE 

CAPSB 

N 

Packed  Stripper  Installed  Capital  Costs,  BENZENE 

CARR 

N 

Total  Pounds  Carbon  Needed,  Time  Base 

ecu 

N 

Carbon  Cycle  Use  BENZENE,  pounds/cycle 

CFPD 

N 

Carbon  Cubic  Feet  Required,  Time  Base 

CMH 

N 

Contaminant  Mole  Fraction,  In 

CMFO 

N 

Contaminant  Mole  Fraction,  Out 

CN 

N 

Column  Number  Lookup  Table  9 

CPD 

N 

Adjusted  Carbon  Needed,  pounds.  Time  Base 

CPDl 

N 

CPD,  Individual  for  Benzene 

CPD2 

N 

CPD,  Individual  for  Trichloroethylene 
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TABLE  C-1.  NOTATION  TABLE  -  DESIGN  AND  COST  ESTIMATING  (CONTINUED) 


NAMES 

n>i 

ITEM  DESCRIPTION  [c] 

CPS 

N 

Carbon  Unit  Design  Size  BENZENE,  pounds 

CPV 

N 

Collision  Function,  Vapor 

CRA 

N 

Carbon  Recycled  Annually  BENZENE,  pounds 

CUP 

N 

Carbon  Units  Purchased  [Number] 

DIB 

T 

Table  DIB,  Area  Name,  BENZENE 

DIT 

T 

Table  DIT,  Area  Name,  TCE 

D2B 

T 

Table  D2B,  Area  Name,  BENZENE  Cost  Estimate 

D2T 

T 

Table  D2T,  Area  Name,  TCE  Cost  Estimate 

D3B 

T 

Table  D3B,  Area  Name,  BENZENE  Cost  Estimate 

D3T 

T 

Table  D3T,  Area  Name,  TCE  Cost  Estimate 

D4B 

T 

Table  D4B,  Area  Name,  BENZENE  Cost  Estimate 

D4T 

T 

Table  D4T,  Area  Name,  TCE  Cost  Estimate 

D5B 

T 

Table  D5B,  Area  Name,  BENZENE  Cost  Estimate 

D5T 

T 

Table  D5T,  Area  Name,  TCE  Cost  Estimate 

D6B 

T 

Table  D6B,  Area  Name,  BENZENE  Cost  Estimate 

D6T 

T 

Table  D6T,  Area  Name,  TCE  Cost  Estimate 

D7B 

T 

Table  D7B,  Area  Name,  BENZENE  Cost  Estimate 

D7T 

T 

Table  D7T,  Area  Name,  TCE  Cost  Estimate 

DAK 

N 

Diffusivity,  m^/second 

DAKA 

N 

Diffusivity,  ft^/second 

DALR 

N 

Derived  Air  Loading  Rate,  (pounds/minute )/ft^ 

DIA 

N 

Stripper  Diameter,  feet 

DIAl 

N 

Stripper  Diameter,  Table  DIB 

DIN 

N 

Stripper  Diameter,  Inches 

DORD 

N 

Derived  Ordinate 

TABLE  C-1.  NOTATION  TABLE  -  DESIGN  AND  CXDST  ESTIMATING  (CONTINUED) 


NAMES 

[bj 

ITEM  DESCRIPTION  [c] 

DORDB 

N 

DORD,  Individual  for  Benzene 

DORDT 

N 

DORD,  Individual  for  TCE 

EXPCATB 

N 

Catalysis  Unit  Annual  Expense,  BENZENE 

EXPCB 

N 

Carbon  Unit  Annual  Expense,  BENZENE 

EXPRB 

N 

Centrifugal  Contactor  Annual  Expense,  BENZENE 

EXPSB 

N 

Packed  Stripper  Annual  Expense,  BENZENE 

FPC 

N 

Flood  Point  Coefficient 

FUELIB 

N 

Carbon  Air  Heater  Fuel,  10^  BTU/year 

FUEL2B 

N 

Catalysis  Unit  Fuel,  10^  BTU^ear 

GCH 

N 

Gravity  Acceleration,  Time  =  Hours 

GCM 

N 

Gravity  Acceleration,  Time  =  Minutes 

GCS 

N 

Gravity  Acceleration,  Time  =  Seconds 

H 

N 

Henry’s  Law  Cbnstants 

HA 

N 

HTU,  Air  Feet 

HAl 

N 

Height  of  Transfer  Unit,  Table  DIB 

HF 

N 

Total  Friction  Loss  Factor 

HPl 

N 

Air  Blower,  Total  Horsepower 

HP2 

N 

Water  Pump,  Total  Horsepower 

HT 

N 

HTU,  Total  Feet 

HTOT 

N 

Total  Height  of  Stripper  Packing,  Feet 

HTOTl 

N 

Stripper  Total  Height,  Table  DIB 

HTR 

N 

Total  Stripper  Height,  Feet 

HW 

N 

HTU,  Water  Feet 

IN 

N 

Table  Index  Column,  Table  8 

INPD 

N 

Total  Pressure  Drop,  inches  of  water 
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TABLE  C-1.  NOTATION  TABLE  -  DESIGN  AND  CX)ST  ESTIMATING  (CXDNTINUED) 


NAMES 

P>1 

ITEM  DESCRIPTION  [c] 

INPDl 

N 

Total  Pressure  Drop,  inches.  Table  DIB 

INPD2 

N 

INPD,  Individual  for  TCE 

INT 

N 

Flooding  Ordinate,  Intermediate  Term  Calculation 

INT2 

N 

Intermediate  Calculation  Factor,  Table  5 

INT3 

N 

Intermediate  Calculation  Factor,  Table  5 

INT4 

N 

Intermediate  Calculation  Factor,  Table  S 

K 

N 

Operating  Ordinate 

K1 

N 

K  Individual  for  Benzene 

K2 

N 

K  Individual  for  Trichloroethylene 

KWl 

N 

Air  Blower,  Total  Kilowatts 

KW2 

N 

Water  Pump,  Total  Kilowatts 

KWS 

N 

Total  Stripper  Electric  Power,  Kilowatts 

KWSl 

•  N 

KWS,  Individual  for  Benzene 

KWS2 

N 

KWS,  Individual  for  TCE 

LPP 

N 

Liquid  Packing  Parameter,  B&F 

LPPl 

N 

Liquid  Packing  Parameter,  B&F 

LPP2 

N 

Liquid  Packing  Parameter,  B&F 

LPP3 

N 

Liquid  Packing  Parameter,  B&F 

LPP4 

N 

Liquid  Packing  Parameter,  B&F 

MASSR 

N 

Mass  Ratio,  Water/Air 

MOLER 

N 

Mole  Ratio,  Water/Air 

MRAT 

N 

Air  Stream  Concentration,  ppm  by  moles 

MRATl 

N 

MRAT,  Individual  for  Benzene 

MRAT2 

N 

MRAT,  Individual  for  TCE 

MSEP 

N 

Molecular  Separations,  Air-Contaminant 
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TABLE  C-1.  NOTATION  TABLE  -  DESIGN  AND  COST  ESTIMATING  (CONTINUED) 


NAMES 

lb] 

ITEM  DESCRIPTION  [c] 

MWA 

N 

Molecular  Weight,  Air 

MWK 

N 

Molecular  Weight,  Contaminant 

MWW 

N 

Molecular  Weight,  Water 

MTU 

N 

Number  of  Transfer  Units 

NTUl 

N 

Number  of  Transfer  Units,  Table  DIB 

ORD 

N 

Flooding  Ordinate 

PAGB 

N 

Pressure  Gradient  for  Benzene,  pounds/^n^  foot 

PAGBI 

N 

Pressure  Gradient  for  Benzene,  inches  water/foot 

PAGT 

N 

Pressure  Gradient  for  TCE,  pounds^n^  foot 

PAGTI 

N 

Pressure  Gradient  for  TCE,  inches  water/foot 

PB 

N 

Pressure  Gradient  BENZENE,  pounds/ft^  foot 

PC 

N 

Plastic  Pall  Rings,  $/ft^ 

PD 

N 

Water  Pipe  Diameter,  feet 

PDIN 

N 

Water  Pipe  Internal  Diameter,  inches 

PGT 

N 

Pressure  Gradient  Data  Array,  Table  8 

PPBD 

N 

Contaminant  Concentration  Change,  ppb 

PPD 

N 

Total  Pressure  Drop,  pounds/in^ 

PT 

N 

Pressure  Gradient  TCE,  pounds/ft^  foot 

PVOL 

N 

Packing  Volume,  ft’ 

PVOLl 

N 

Packing  Volume,  Table  DIB 

RCP 

N 

Computation  Term,  Vapor  HTU 

RN 

N 

Row  Number  Lookup  Table  9 

RRPMM 

N 

Contaminant  Removal  Rate,  pound  moles/minute 

RRPPM 

N 

Contaminant  Removal  Rate,  pounds/minute 

SELB 

N 

Benzene  Value  Finder,  Table  8 
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TABLE  C-1.  NOTATION  TABLE  -  DESIGN  AND  CX)ST  ESTIMATING  (CONTINUED) 


ITEM  DESCRIPTION  [c] 


Ordinate  Selector  BENZENE,  Table  9 


Ordinate  Selector  TCE,  Table  9 


Trichloroethylene  Value  Finder,  Table  8 


Control  Setup  for  Sensitivity  Analysis 


Schmidt  Number,  Vapor 


Schmidt  Number,  Weight 


Stripper  Weight,  Table  DIB 


Treybal  Abscissa,  Selected  for  Table  8 


Table  1,  Area  Name,  Design 


Table  2,  Area  Name,  Design 


Table  3,  Area  Name,  Design 


Table  4,  Area  Name,  Design 


Table  5,  Area  Name,  Design 


Table  6,  Area  Name,  Design 


Table  7,  Area  Name,  Design 


Table  8,  Area  Name^  Design 


Table  9,  Area  Name,  Design 


Table  10,  Area  Name,  Design 


Table  11,  Area  Name,  Design 


Table  12,  Area  Name,  Design 


Celsius  Temperature 


Total  Pounds  Contaminant  Stripped,  Time  Base 


Kelvin  Temperature 


Treybal  Ordinate  Table  9 


Total  Pipe  Length  for  Water,  feet 
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TABLE  C- J  NOTATION  TABLE  -  DESIGN  AND  COST  ESTIMATING  (CONCLUDED) 


NAMES 

[b] 

ITEM  DESCRIPTION  [c] 

TPW 

N 

Total  Pounds  Water,  Time  Base 

TR 

N 

Rankine  Temperature 

VPP 

N 

Vapor  Packing  Parameter,  B&F 

VPPl 

N 

Vapor  Packing  Parameter,  B&F 

VPP2 

N 

Vapor  Packing  Parameter,  B&F 

VPP3 

N 

Vapor  Packing  Parameter,  B&F 

VPP4 

N 

Vapor  Packing  Parameter,  B&F 

W1 

N 

Air  Blower,  Total  Watts 

W2 

N 

Water  Pump,  Total  Watts 

WAK 

N 

Water  Diffusivity,  m^/second 

WAKA 

N 

Water  Diffusivity,  ft^/second 

WCFM 

N 

Water  Feed  Rate,  ftVminute 

WD 

N 

Water  Density,  pounds/ft^ 

WLR 

N 

Water  Loading  Rate,  (pounds/ft^)/hour 

WPM 

N 

Water  Feed  Rate,  pounds/minute 

WPMM 

N 

Water  Feed  Rate,  pound  moles/minute 

WRAT 

N 

Air  Stream  Concentration,  ppm  by  weight 

WRATl 

N 

WRAT,  Individual  for  Benzene 

WRAT2 

N 

WRAT,  Individual  for  TCE 

WVl 

N 

Water  Viscosity,  Centipoise 

WV2 

N 

Water  Viscosity,  pounds/foot  second 

WVSI 

N 

Water  Viscosity,  Pa  second 

Z 

N 

Total  Vertical  Lift  for  Water,  feet 

[a]  +  =  Designation  of  a  variable  design  or  cost  input  parameter. 

[b]  A  =  Alpha  entry,  N  =  Numeric  entry,  T  =  Table  name  for  printing. 

[c]  Descriptive  Comments.  Names  used  in  Tables  D-DIT  through  D-D7T  for 
TCE  are  not  included  since  they  are  parallel  to  those  used  in  BENZENE 
Tables  D-D  IB  through  D-D7B. 
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APPENDIX  D 


DESIGN  AND  COST-ESTIMATING  SPREADSHEET 

VOC-2.100 

The  following  26  tables  make  up  spreadsheet  VOC-2.1(X),  which  is  written  in  MULTIPLAN 
2.01.  The  file  may  be  exported  upward  into  MP  3.0  or  MP  4.0  without  difficulty.  Tables  D-1 
through  D-12  deal  with  the  design  parameters  of  a  packed  column,  a  carbon  adsorption  off-gas 
treatment  system,  and  a  catalytic  oxidation  off-gas  treatment  system.  These  12  tables 
simultaneously  handle  TCE  and  Benzene,  the  latter  used  as  a  marker  for  jet  fuel.  The  seven 
tables  D-DIB  through  D-D7B  are  sensitivity  analysis  output  summaries  for  the  Benzene  case,  and 
the  seven  tables  D-DIT  through  D-D7T  are  the  same  for  the  TCE  case.  Tables  D-DIB,  D-D2B, 
D-D  IT,  and  D-D2T  summarize  requisite  design  information  for  use  in  single-variant  sensitivity 
analyses.  Tables  D-D3B  through  D-D5B  and  D-D3T  through  D-D5T  summarize  the  fixed  capital 
and  the  annual  operating  cost  estimates  for  the  packed  column  stripper,  the  carbon  adsorption 
off-gas  treatment,  and  the  catalytic  oxidation  off-gas  treatment.  Tables  D-D6B  and  D-D6T 
consolidate  the  fixed  capital  and  the  annual  operating  cost  estimates  for  the  packed  column 
stripper  options.  Tables  D-D7B  and  D-D7T  consolidate  these  for  the  centrifugal  contactor 
option. 

Fixed  capital  and  annual  operating  (non-capital)  cost  estimates  have  been  adjusted  to  1990 
dollars  in  the  following  tables.  Continuing  inflation  adjustments  for  the  period  of  the  operating 
lifetime  are  handled  in  a  separate  spreadsheet,  VOC-2.300,  which  is  described  in  Appendix  E. 
Further  details  will  be  found  in  the  related  body  of  the  text. 
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DESIGN  INPUT  TABLE 
Indapandant  0«al9n  Variabla  Inputs 
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Indicataa  original  location  of  a  NAMED  INDEPENDENT  VARIABLE  input  (see  Table  1). 
Indicates  original  location  of  a  NAMED  DEPENDENT  VARIABLE. 

Indicates  a  NAMED  VALUE,  transferred  from  an  earlier  location. 
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.100  DIAMETER  CALCULATIONS  FOR  PACKED  STRIPPER 
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[TABS]  TABLE  D-5 

VOC-2.100  VAPOR  HTU  CALCULATIONS  FOR  PACKED  STRIPPER 
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This  diameter  value  is  used  for  calculation  of  the  recjuired  volume  of  packing,  and  for 

the  cross  sectional  area  of  the  packed  column  stripper.  BUT  ...  an  artificial 

dieuneter  of  2  feet  appears  in  the  calculation  of  the  HTU  [HA]  in  this  table.  [Hylton,  ORNL) 
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[TABS]  TABLE  D-8 

VOC-2.100  NUMBER  OF  TRANSFER  UNITS  CALCULATIONS  FOR  PACKED  STRIPPER 
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TOTAL  HEIGHT  OF  PACKING,  FT  (HTOT)  53.67  56.36 
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OPERATING  VARIABLES 
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(TABlll  TABLE  D-11 

VOC-2.100  AIR  BLOWER  &  WATER  PUMP  CALCULATIONS  FOR  PACKED  STRIPPER 
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TOTAL  STRIPPER  POWER,  kW,  (KWSB],  [KWST]  (KWSJ  14.94  14.35 
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APPENDIX  E 


OPERATING  LIFETIME  FINANCIAL  ANALYSIS  SPREADSHEET 


VOC-2.3QO 

The  following  eight  tables  make  up  spreadsheet  VOC-2.300  written  in  MULTIPLAN  2.01.  The 
file  may  be  exported  upward  into  MP  3.0  or  MP  4.0  without  difficulty.  Table  E-1  is  for  inputs  to  the 
construction  phase,  and  Table  E-4  is  for  inputs  to  the  operating  lifetime  phase.  For  the  analyses  in 
this  report,  the  construction  phase  is  assumed  to  be  zero  years,  so  that  Tables  E-1  and  E-4  are 
essentially  identical.  Table  E-8  contains  the  output  information,  in  terms  of  lifetime  current  dollars 
cost  per  1000  gallons  of  groundwater  processed.  Annual  adjustments  for  inflation  are  included,  using 
the  input  inflation  projection  shown  in  Table  E-4.  The  analyses  are  handled  as  100  percent  debt 
financing  to  simulate  a  government-owned  and  operated  activity.  At  100  percent  debt  financing  and 
at  zero  salvage  on  the  capital  equipment,  the  lifetime  processing  costs  per  1000  gallons  of  water 
should  be  identical  for  a  cash  flow  analysis  and  for  a  depreciation  analysis.  This  identity  is  shown  in 
Table  E-8.  For  cash  flow  analysis,  the  costs  include  the  total  capital  costs  and  the  sum  of  the  yearly 
operating  expenses.  For  the  depreciation  analysis,  the  operating  expenses  are  the  same  as  for  the 
cash  flow  analysis.  The  capital  costs,  though,  are  the  sum  of  the  yearly  depreciation  amounts.  Thus, 
if  the  capital  equipment  is  depreciated  to  a  zero  salvage  value,  the  two  analyses  are  identical.  The 
operating  lifetime  financial  analysis  is  discussed  in  the  related  body  of  the  text 

The  lifetime  processing  costs  shown  in  Table  E-8  are  for  a  particular  set  of  input  variables.  In 
Section  V  of  the  text,  exarooles  were  presented  which  showed  how  processing  costs  were  affected  by 
variations  in  certain  variables.  Additional  examples  of  such  sensitivity  analyses  are  presented  in 
Figures  E-1  through  E-4,  which  show  the  effects  of  variations  in  overhead  rate,  carbon  recycle 
interval,  contaminant  removal  efficiency,  and  equipment  installation  factor. 
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(PF],  PF  - >  PRINT  FILE  MACRO,  FOR  ALL  TABLES  IN  SPREADSHEET 


BEFORE  &  AFTER  INFLATION  ADJUSTMENTS, 
WITH  DEBT  FRACTION  ADJUSTMENT  AS  WELL 
BUT  ...  WITH  NO  INTEREST  INCLUDED. 
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[D]  *  1.000  (INOJ  »  0.098 

[DO]  »  1.000  (IFOJ  -  0.048 


DOLLAR  COST  PER  1,000  GALLONS,  BASED  UPON  CASHFLOW  ANALYSIS  »  1.5208 

DOLLAR  COST  PER  1,000  GALLONS,  CONTRIBUTION  OF  NON-CAPITAL  ANNUAL  EXPENSES  -  1.4736 

DOLLAR  COST  PER  1,000  GALLONS,  CONTRIBUTION  OF  FIXED  CAPITAL  DEPRECIATION  -  0.0472 
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Operating  Lifetime  Processing  Costs  as  a  Function  of  Overhead  Rate  for  an  Air  Stripper  System  Featuring 
Traditional  Packed  Tower  for  TCE  Removal  from  Groundwater. 
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Figure  E-2.  Operating  Lifetime  Processing  Costs  as  a  Function  of  Carbon  Recycle  Interval  for  an  Air  Stripping  System  Featuring 
a  Traditional  Packed  Tower  for  TCE  Removal  from  Groundwater. 


Stripper  System  Featuring  a  Traditional  Packed  Tower. 
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